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keying  (BPSK).  Matlab  was  used  to  simulate  the  design,  which  was  tested  in  both  an 
additive  white  Gaussian  noise  (AWGN)  channel  and  in  a  slow  fading  frequency 
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the  maximal  ratio  combiner  (MRC)  receiving  technique  with  perfect  knowledge  of 
channel  state  information  (CSI).  The  theoretical  perfonnance  was  derived  for  both 
channels  and  was  compared  with  the  simulated  results. 
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EXECUTIVE  SUMMARY 


In  modern  military  and  commercial  wireless  communications,  the  demand  for 
high  speed  and  reliable  communication  within  the  constraints  of  limited  radio  frequency 
spectrum  and  power,  are  the  prime  technical  criteria  for  communication  systems.  To 
obtain  a  higher  data  rate  at  an  acceptable  bit  error  rate,  larger  bandwidth  is  required.  To 
mitigate  severe  fading  channel  conditions,  a  higher  transmitted  power  level  is  required. 
Multiple-input  multiple-output  (MIMO)  communication  systems  have  the  potential  to 
provide  increased  capacity  and  reliability  without  increasing  the  bandwidth  or  transmitted 
power.  MIMO  systems  exploit  time  and  spatial  diversities  by  employing  multiple 
antennas  at  the  transmitter  and  receiver.  MIMO  systems  may  help  to  enhance  the 
robustness  of  military  communication  systems  under  the  worst  operating  conditions. 

Multicarrier  communication  in  the  form  of  orthogonal  frequency  division 
multiplexing  (OFDM)  has  been  adopted  in  several  wireless  communications  standards 
due  to  its  ability  to  mitigate  severe  multipath  conditions,  its  bandwidth  efficiency  and  its 
simplicity  of  implementation.  The  combination  of  OFDM  and  MIMO  technologies  holds 
the  promise  of  increasing  future  communications  demands. 

The  main  objective  of  this  thesis  was  to  investigate  the  fundamentals  of  MIMO 
systems  with  a  multicarrier  delay  diversity  modulation  (MDDM)  technique.  A  simple 
model  was  designed  to  incorporate  MDDM  in  multiple-input  single-output  (MISO)  and 
MIMO  systems.  This  design  was  simulated  and  analyzed  to  demonstrate  its  perfonnance. 
The  system  was  implemented  with  binary  phase  shift  keying  (BPSK)  in  Matlab  and  was 
tested  in  both  an  additive  white  Gaussian  noise  (AWGN)  channel  with  no  fading  and  a 
slow  multipath  fading  channel  with  AWGN.  The  receiver  design  was  based  on  the 
maximal  ratio  combining  (MRC)  technique  with  the  assumption  of  perfect  knowledge  of 
channel  state  information  (CSI)  at  the  receiver  end.  The  simulated  performance  results 
and  theoretical  analysis  results  were  compared  with  the  conventional  single-input  single¬ 
output  (SISO)  system  results.  The  perfonnance  metric  of  bit  error  probability  versus 


Eb  /  N0  (energy  per  bit  to  noise  power  spectral  density  ratio)  was  used.  To  establish  a 
fair  comparison,  the  transmitted  power  for  the  SISO,  MISO  and  MIMO  systems  was 
maintained  equal. 

The  results  showed  that  the  designed  MISO  and  MIMO  system  performed  within 
expected  parameters  of  the  theoretical  analysis  in  both  the  AWGN  channel  with  no 
fading  and  the  multipath  fading  channel  with  AWGN.  The  comparison  of  perfonnances 
in  the  AWGN  channel  with  no  fading  showed  that  the  MISO  system  performed  better 
than  the  SISO  system  for  low  Eb  /  N0  values  up  to  6.5  dB  and  the  perfonnance  of  the 
MISO  system  was  poorer  for  higher  Eb  /  N0  values.  The  performances  of  the  MIMO 
systems  were  better  than  that  of  the  SISO  system  for  ah  values  of  Eb  /  N0  and  all 
systems  studies  herein.  The  MIMO  systems  with  two  receive  antennas  and  three  receive 
antennas  outperformed  a  SISO  system  by  1.7  dB  and  3.4  dB  less  transmit  power  required 
respectively  for  equal  performance.  For  the  multipath  fading  channel  with  AWGN,  the 
MISO  and  MIMO  systems  were  able  to  achieve  significant  advantage  over  a  SISO 
system. 


I.  INTRODUCTION 


A.  BACKGROUND 

One  of  the  major  challenges  facing  modern  communications  is  to  satisfy  the  ever 
increasing  demand  of  high  speed  reliable  communications  with  the  constraints  of 
extremely  limited  frequency  spectrum  and  limited  power.  Wireless  communications 
systems  like  cellular  mobile  communications,  internet  and  multimedia  services  require 
very  high  capacity  to  fulfill  the  demand  of  high  data  rates.  These  systems  must  achieve 
the  desired  reliability  within  the  limits  of  power  and  frequency  spectrum  availability, 
often  in  severe  channel  environments.  They  need  to  overcome  signal  scattering  and 
multipath  effects,  especially  in  densely  populated  urban  areas.  For  many  military 
commutation  systems,  reliable  communication  is  to  be  achieved  with  low  probability  of 
detection  and  interception  even  in  hostile  jamming  environments. 

The  solutions  to  achieve  high  capacity  with  reliability  could  include  time, 
frequency  and  space  diversity.  Wireless  communication  systems  with  multiple  transmit 
and  multiple  receive  antennas  can  provide  high  capacity  at  low  probability  of  bit  error 
with  extremely  low  power,  even  in  dense  scattering  and  multipath  environments.  These 
multiple-input  multiple-output  (MIMO)  systems  with  appropriate  space-time  codes  have 
been  an  area  of  recent  research  as  they  hold  the  promise  of  ever  increasing  data  rates.  The 
capacity  of  a  MIMO  system  can  be  increased  linearly  by  increasing  the  number  of 
transmit  and  receive  antennas.  The  applications  of  MIMO  systems  in  a  frequency 
selective  channel  require  equalization  and  other  techniques  to  compensate  for  frequency 
selectivity  of  the  channel,  which  add  to  the  complexity  of  these  systems.  In  recent  years, 
orthogonal  frequency  division  multiplexing  (OFDM)  has  been  widely  used  in 
communications  systems  to  operate  in  frequency  selective  channels  including  several 
wireless  communication  standards.  Communication  systems  with  a  MIMO-OFDM 
combination  can  significantly  improve  capacity  and  reliability  by  exploiting  the 
robustness  of  OFDM  to  fading,  enhanced  by  adding  more  diversity  gain  via  space  time 
codes.  In  this  research,  delay  diversity  in  combination  with  a  MIMO-OFDM  system  is 
the  primary  focus  and  is  referred  to  as  multicarrier  delay  diversity  modulation  (MDDM). 
[1,2] 
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B.  OBJECTIVE  AND  METHODOLOGY 

The  main  objective  of  this  research  was  to  investigate  MIMO-OFDM  systems  by 
using  a  cyclic  delay  diversity  technique  as  a  space  time  code  and  compare  its 
performance  to  conventional  single-input  single-output  (SISO)  systems.  The  first  step  to 
achieve  the  objective  was  to  study  the  fundamentals  of  OFDM  systems.  Then,  the  OFDM 
technique  was  applied  to  a  multiple-input  single-output  (MISO)  system  and  MIMO 
systems  with  cyclic  delay  diversity.  A  MISO  system  and  two  MIMO  systems,  using  a 
MDDM  scheme  with  two  transmit  antennas,  were  simulated  in  Matlab  with  one,  two  and 
three  receive  antennas,  respectively.  Multicarrier  delay  diversity  was  implemented  with 
binary  phase  shift  keying  (BPSK)  and  was  realized  in  equivalent  baseband  form  to 
facilitate  comparison  with  published  theory.  The  design  was  simulated  both  with  an 
additive  white  Gaussian  noise  (AWGN)  channel  with  no  fading  and  then  with  a  multipath 
faded  channel  with  AWGN.  Theoretical  perfonnance  and  simulated  system  performance 
were  compared  with  reference  to  the  single-input  single-output  (SISO)  system. 

C.  RELATED  RESEARCH 

MIMO  systems  promise  high  efficiency  in  providing  low  probability  of  bit  error 
without  increasing  the  transmitted  power  or  the  bandwidth  [1].  Therefore,  research  in  this 
area  has  been  very  active  during  recent  years.  The  perfonnance  of  MIMO  systems 
applying  techniques  such  as  delay  diversity,  space-time  block  codes  (STBC)  and  space- 
time  trellis  codes  (STTC)  has  led  to  the  development  of  several  practical  systems  [1,  2]. 
Numerous  studies  have  been  performed  to  realize  and  investigate  the  performance  of 
these  systems.  Space  time  codes  were  originally  developed  for  flat  fading  channels. 
Later,  incorporation  of  OFDM  in  these  systems  has  provided  a  solution  to 
communications  in  the  frequency  selective  channel  [1].  Delay  diversity  was  the  first 
diversity  technique  introduced  for  MIMO  systems  and  was  described  in  [3]  and  [4].  The 
use  of  delay  diversity  with  OFDM  was  proposed  in  [5]  for  a  flat  fading  channel  and  a 
cyclic  delay  diversity  approach  with  OFDM  was  recommended  for  the  frequency 
selective  fading  channel  in  [6]  which  has  been  further  investigated  in  [2]. 

D.  THESIS  ORGANIZATION 

This  thesis  is  organized  in  five  chapters.  Chapter  II  introduces  the  fundamentals 
of  MIMO  systems  and  the  multicarrier  delay  diversity  modulation  scheme.  It  discusses 
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the  capacity  of  MIMO  systems,  space  time  coding  and  use  of  time  delay  diversity  as  a 
space  time  coding  technique  in  MIMO  systems.  Chapter  III  describes  the  modeling  and 
simulation  of  the  MDDM  transmitter  and  receiver.  Chapter  IV  analyzes  the  performance 
of  this  modulation  technique  both  in  an  AWGN  channel  with  no  fading  and  in  a  slow 
fading  multipath  channel  with  AWGN.  Chapter  V  reviews  the  summary  of  the  work  done 
with  results  and  includes  recommendations  for  future  study.  Matlab  code  is  attached  as  an 
Appendix. 
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II.  MIMO  SYSTEMS  AND  MULTICARRIER  DELAY  DIVERSITY 


The  objective  of  this  chapter  is  to  establish  the  basic  understanding  of  MIMO 
systems,  space  time  coding  and  the  application  of  MDDM  in  MIMO  systems.  An 
overview  of  the  MIMO  model  is  presented.  The  fundamentals  of  OFDM  and  its 
implementations  with  the  discrete  Fourier  transform  (DFT)  are  discussed.  Then,  the 
implementation  of  cyclic  delay  diversity  with  OFDM  within  the  context  of  MIMO 
systems  is  presented.  Finally,  the  multipath  flat  fading  channel  and  the  maximum  ratio 
combining  receiver  are  discussed. 

A.  MULTIPLE-INPUT  MULTIPLE-OUTPUT  (MIMO)  SYSTEMS 

In  order  to  facilitate  the  understanding  of  multiple-input  multiple-output  systems, 
single-input  single-output  (SISO)  systems,  single-input  multiple-output  (SIMO)  systems 
and  multiple-input  single-output  (MISO)  systems  models  are  discussed  briefly.  In  this 
thesis,  all  signals  and  models  are  represented  in  complex  baseband  equivalent  form  to 
facilitate  analysis. 

1.  Single-Input  Single-Output  System 

The  SISO  system  model  is  shown  in  Figure  1.  The  signal  transmitted  from  the 
antenna  is  denoted  as  x(t) .  The  signal  received  at  the  receiving  end,  r(t),  passes  through 
the  channel  with  impulse  response  h(t )  in  an  additive  white  Gaussian  noise  (AWGN) 
environment.  It  is  assumed  that  the  bandwidth  of  the  signal  is  small  enough  such  that  the 
frequency  response  of  the  channel  is  flat  and  the  channel  response  can  be  given  as 

h  =  \h\e4  (2.1) 

where  i  =  VM  .  The  relationship  between  the  transmitted  signal  and  receive  signal  is 
given  by 

r(V)  =  J  x{x)h(t-x)di  +  n(t) 

=  |  x{r^\h\e1<pS{t-T)dT +  n{t)  (2.2) 

=  x{t)h  +  n(t) 
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where  S(t)  is  the  Dirac  delta  function  and  n(t)  is  the  AWGN.  The  received  signal  is  the 
transmitted  signal  convolved  with  the  channel  impulse  response  plus  added  noise.  [7] 


Transmitter 


Figure  1  SISO  System  (After  Ref.  [7]) 

2.  Single-Input  Multiple-Output  System 

A  single-input  multiple-output  (SIMO)  system  with  one  transmit  antenna  and 
multiple  (J)  receive  antennas  is  illustrated  in  Figure  2. 


Figure  2  SIMO  System  (After  Ref.  [7]) 

The  received  signals  at  the  receiver  end  can  be  represented  as  the  set  of  linear  equations 


rx(f)  =  hl{t)x(t)  +  n{t) 
r2(t )  =  h2(t)x{t)  +  n2(t) 

(2.3) 

rJ(t )  =  hJ(t)x(t)  +  nJ(t ) 
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where  rJ(t),  hJ(t)  and  nJ(t),  j  =  1,2, 3,...,  J  represent  the  received  signal,  channel 
impulse  response  and  noise,  respectively,  at  the  j"th  receive  antenna.  If 

x(t),  h(t )  and  r(t)  are  sampled  at  the  rate  of  one  sample  per  symbol,  then  they  can  be 
represented  as  x,  h,  r  .  The  received  signal  can  also  be  represented  in  form  of  vectors. 

x  =  [x'  x2  •••  for  x1  =  x2  = ...  =  xJ  (2.4) 

h  =  ^/z1  h2  •••  hJ  J 

n  =  [V  n2  •••  nJ  J  (2.5) 

Now,  the  received  vector  is  represented  as 

r  =  h.x  +  n  (2.6) 

where  x,  h  and  n  are  transmission,  channel  and  noise  vectors,  respectively  and  operator 
' . '  denotes  element  by  element  multiplication. 

3.  Multiple-Input  Single-Output  System 

A  multiple-input  single-output  (MISO)  system  with  multiple  (L)  transmit 

antennas  and  one  receive  antenna  is  illustrated  in  Figure  3.  The  receive  antenna  receives  a 
sum  of  all  the  signals  transmitted  by  each  antenna  and  can  be  represented  as 

r  =  hlxl  +  h2x2  +  ...  +  hLxL  +  n  (2.7) 

where  xl  and  hl,  I  =  1.2, 3, are  the  transmitted  signal  and  channel  response  from 
transmit  antenna  l  to  the  receive  antenna  and  n  is  the  AWGN.  Equation  (2.7)  can  also  be 
written  as 

L 

r  =  ^  hlxl  +  n .  (2.8) 

;=i 
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Figure  3  MISO  System  (After  Ref.  [7]) 

4.  Multiple-Input  Multiple-Output  System 

A  multiple-input  multiple-output  (MIMO)  system  with  multiple  (L)  transmit 
antennas  and  multiple  (J)  receive  antennas  is  illustrated  in  Figure  4.  MIMO  systems  are 
the  focus  of  this  thesis.  Therefore,  the  MIMO  model  is  discussed  in  detail  in  this  section. 
The  representation  of  the  model  is  largely  based  on  [1,  8].  For  a  faded  channel,  it  is 
assumed  that  channel  responses  from  each  transmit  antenna  to  each  receive  antenna  are 
independent. 


Figure  4  MIMO  System  (After  Ref.  [7]) 

The  received  signal  at  any  receive  antenna  is  the  sum  of  all  the  signals  transmitted  by 
each  transmit  antenna  passing  through  the  respective  channel  and  the  AWGN. 
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rj  =  h'Jx'  +  h2Jx2  + . . .  +  hLjxL  +  ri 


:  ~h'J  h2J  ■■■hLj~ 


+  n 


(2.9) 


li!1  is  the  channel  response  from  transmit  antenna  l  to  receive  antenna  j  where 
l  =  1,2,3 and  j  =  1,2, 3,...,  J .  The  system  therefore  can  be  represented  in  matrix 
form  as 


1 

r 

~hn 

h21  ■  ■ 

■hn  " 

i 

X 

1 

n 

2 

r 

h21 

h22  ■ 

•• hL 2 

x2 

n1 

(2.10) 

v'J 

hji 

hJ1- 

■■hJL 

XL 

+ 

nJ 

r  =  Hx  +  n  (2.11) 

where  r  and  n  represent  the  J  x  1  received  signal  and  noise  column  vectors,  H  is  a 
J  x  L  complex  channel  matrix  and  x  is  the  Lxl  transmitted  column  matrix.  [1] 

All  the  elements  in  the  channel  matrix  are  considered  independent  identically 
distributed  (IID)  complex  Gaussian  random  variables  (GRVs),  and  similarly,  the 
elements  of  the  noise  vector  are  also  complex  Gaussian  random  variables.  According  to 
information  theory,  the  optimum  distributions  for  transmitted  signals  are  also  Gaussian. 
Therefore,  it  is  considered  that  all  elements  of  x  are  zero  mean  complex  GRVs. 
Therefore,  the  covariance  matrix  of  the  transmitted  signals  is  the  same  as  the 
autocorrelation  matrix  and  it  is  a  diagonal  matrix  [1,8] 


(2.12) 


where  x  r  denotes  the  transpose  and  component- wise  complex  conjugate  of  the 
transmitted  matrix  x  .  [1] 
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The  total  transmitted  power,  P ,  is  the  sum  of  all  the  diagonal  elements  of  the 
autocorrelation  matrix.  To  facilitate  the  analysis  for  the  MIMO  system,  assume  that  all 
the  transmit  antennas  transmit  equal  power  [8], 


P  =  E 


*(R») 


(2.13) 


P‘  =  E 


P 

~L 


(2.14) 


where  P1  is  the  average  power  transmitted  from  antenna  /.  An  AWGN  channel  is 
considered,  and  according  to  information  theory,  the  optimum  distribution  for  the 
transmitted  signal  is  also  Gaussian  [1],  Therefore,  the  elements  of  x  as  stated  in  [1]  are 
also  considered  independent  and  identically  distributed  (IID)  Gaussian  variables  with 
zero  mean.  Then,  the  autocorrelation  matrix  can  be  written  as 

R„=fT  (215) 


where  I,  is  the  identity  matrix  of  size  Lx  L . 

It  is  further  assumed  that  the  channel  matrix  is  fixed  at  least  for  the  duration  of 
one  symbol  period  and  there  is  no  attenuation  due  to  path  loss  and  no  amplification  due  to 
antenna  gain.  In  other  words,  each  receive  antenna  receives  the  total  transmitted  power 
regardless  of  its  branch.  Thus,  the  nonnalization  constraint  for  the  elements  of  channel 
matrix  H  for  fixed  coefficient  can  be  represented  as  [1] 

L  2 

Y\hi'\=L  for  j  =  \,2,...J .  (2.16) 

i=\ 

For  a  faded  channel,  the  channel  matrix  elements  are  random  variables  and  the 
nonnalization  constraint  will  apply  to  the  expected  value  of  Equation  (2.16).  This 
nonnalization  constraint  is  required  for  a  fair  comparison  with  SISO  systems  with  equal 
power  transmitted  [8].  It  is  also  assumed  that  the  channel  impulse  response  at  that  time, 
referred  to  as  the  channel  state  information  (CSI),  is  perfectly  known  at  the  receiver  by 
sending  training  symbols.  [1] 
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The  elements  of  the  noise  vector  n  are  considered  IID  complex  Gaussian  random 
variables  with  zero  mean  and  the  variance  of  <j02  /  2  for  both  real  and  imaginary  parts. 
Independence  of  the  noise  elements  and  zero  mean  imply  that  the  autocorrelation  and 
covariance  matrix  are  the  same  diagonal  matrix.  Therefore,  the  covariance  matrix  Cnn  of 
the  noise  vector  can  be  represented  as  [1,  8] 

C„.-R.„=E[nn-r]  (2.17) 

R„=<*X  (2-18) 


The  average  signal  power  at  receive  antenna  j  with  assumed  fixed  channel  coefficients 
can  be  given  by  [1,  8] 


pi  _  p  ri  ri* 


=  E 


=  E 


<  mr  m 

*  X 


YJhl,xlYJh'"r 

1=1  m= 1 

YjYjh'Jx'h 

1=1  m= 1 

=  YJYJh,ih,rE[xixm *] 

/=1  m= 1 

=±i*°r  e 


(2.19) 


where  (  )*  denotes  complex  conjugate.  Substitution  of  Equations  (2.15)  and  (2.16)  into 
Equation  (2.19)  yields 


(2.20) 


Then,  the  average  signal-to-noise  ratio  (SNR)  at  each  receive  antenna,  represented  by  y , 
is  given  by  [1,  8] 


r  = 


(2-21) 


Similarly,  the  autocorrelation  matrix  for  the  received  signal  can  be  represented  as  [8] 
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Rrr=E[rr*r\ 


(2.22) 


R„  =E[(Hx  +  n)(Hx  +  n)‘r 
=  Er(Hx  +  n)((Hxf  +nT) 

1  V  /J  (2.23) 

=  E[Hx(Hx)*r  +Hxn*r  +n(Hxf  +nn*r 

=  E  Hx(Hx)  r  +  E [Hxn*r ]  +  E  (Hxfn  +E[nn‘r], 

By  using  the  identity  of  transposition  of  a  product  of  matrices  [9]  as  follows 

(Hx)*r=x*rH*r.  (2.24) 

Equation  (2.23)  can  be  written  as 

R)r  =E^Hxx*rH*r^  +  E^Hxntr^  +  E^xtrH‘’rn^  +  E^nn~r^] .  (2.25) 


If  it  is  assumed  that  the  channel  coefficients  are  detenninistic  and  signal  matrix  x  and 
noise  matrix  n  are  independent  with  zero  mean.  Equation  (2.25)  yields 


R„.  =  HE xx  rjH  r+E[nn  r 


=  HR  II  +  R  . 


(2.26) 


B.  SPACE  TIME  CODING 


Space  time  coding  is  a  technique  to  achieve  higher  diversity  at  the  receiver  end  to 
mitigate  multipath  fading  without  increasing  the  transmitted  power  or  bandwidth.  Space 
time  coding  holds  the  promise  to  maximize  the  system  capacity.  The  system  capacity  is 
defined  in  [1],  “The  maximum  possible  transmission  rate  such  that  the  probability  of 
error  is  arbitrarily  small.”  The  capacity  of  SISO  system  is  given  by  Shannon’s  capacity 
equation  [10] 


C  =  W\og2{l  +  SNR) 

(  P  \ 
C  -W  log2  1  +  — 

V  cr  ) 


(2.27) 


where  C  ,  W ,  P  and  a2  represent  capacity,  bandwidth,  average  signal  power  and  average 

noise  power,  respectively.  The  capacity  of  a  MIMO  system  in  a  flat  fading  channel  with 

perfect  channel  state  infonnation  is  given  by  [1] 
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(2.28) 


C  =  W  log2 


(  P\ 

1  +  LJ  — - 

v  °  ) 


Space  time  coding  techniques  designed  appropriately  with  MIMO  systems  have 
the  potential  to  achieve  the  channel  capacity  in  Equation  (2.28).  Space  time  coding 
provides  the  diversity  both  in  time  and  space  to  achieve  the  higher  performance  at 
reduced  transmitted  power  and  without  bandwidth  expansion.  Space  time  coding 
techniques  can  be  classified  into  two  main  categories.  The  first  category  provides  power 
efficiency  without  compromising  the  performance  such  as  delay  diversity,  space-time 
block  codes  (STBC)  and  space-time  turbo  trellis  codes  (STTC).  The  other  category,  such 
as  Bell  Labs  layered  space-time  technology  (BLAST),  increases  the  data  rates  with  the 
use  of  bandwidth  efficient  modulation  schemes  [2]. 

The  performance  of  a  MIMO  system  can  be  further  improved  by  applying 
forward  error  correction  coding  (FEC)  with  optimum  interleaving  at  the  cost  of  reduced 
data  rate  or  increased  bandwidth.  FEC  coding  gain  can  be  achieved  without  sacrificing 
the  data  rate  or  bandwidth  by  designing  space  time  coding  technique  with  higher  rate 
modulation  schemes.  [1] 

This  thesis  is  focused  on  achieving  multicarrier  delay  diversity  gain.  The 
incorporation  of  error  control  coding  and  interleaving  is  left  for  future  work. 

C.  MULTICARRIER  DELAY  DIVERSITY  IN  MIMO  SYSTEMS 

The  delay  diversity  technique  was  the  first  approach  proposed  for  MIMO  systems 
[2].  In  this  scheme,  delayed  versions  of  the  same  signal  are  transmitted  by  multiple 
antennas.  This  simple  delay  diversity  was  originally  suggested  for  flat  fading  channels 
[11].  This  scheme  has  the  inherent  problem  of  increasing  frequency  selectivity  caused  by 
the  delay  diversity.  Full  diversity  cannot  be  achieved  without  equalization  [12]  and 
equalization  for  MIMO  systems  is  very  difficult  due  to  the  large  number  of  channels. 
Thus,  the  receiver  design  becomes  much  more  complicated.  Orthogonal  frequency 
division  multiplexing  (OFDM)  with  delay  diversity  is  another  approach  to  make  good  use 
of  frequency  selectivity  of  the  delay  diversity.  The  OFDM  scheme  with  delay  diversity 
has  a  limitation  as  an  increase  in  the  number  of  transmit  antennas  requires  an  increase  in 
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the  guard  interval  at  the  expense  of  bandwidth.  If  the  guard  interval  is  not  as  large  as  or 
larger  than  the  delay  spread  of  the  channels,  then  it  will  cause  inter-symbol  interference 
(ISI)  [12]. 

The  OFDM  scheme  can  be  easily  implemented  by  using  the  inverse  discrete 
Fourier  transform  (IDFT),  which  can  be  efficiently  computed  by  the  inverse  fast  Fourier 
transfonn  (IFFT).  The  same  information  can  be  translated  back  by  the  DFT  operation. 
DFT  and  IDFT  as  defined  in  [13]  are 

N- 1 

DFT  {£[???.]}  =X[k]  =  £a;[m]  e-^kn/N  for  k  =  0,1,2,...,#  -  1  (2.29) 

n=0 

IDFT {#[&]}  =  xl ml  =  —  ^X[k]  ^nkn/N  for  rn  =  0,1,2,...,#  -  1 .  (2.30) 

X  k=o 

Discrete  Fourier  transforms  have  a  property  that  any  circular  shift  in  the  time  domain 
results  in  a  phase  shift  in  the  frequency  domain  [13] 

DFT [x[(m-D)N]}  =  e-t2!cW/NX[k]  for  k  =  0,1,2,...,#  - 1  (2.31) 

where  D  denotes  the  delay  in  the  time  index  and  (n  )  v  denotes  n  modulo  #  .  This  cyclic 
delay  property  of  the  Fourier  transforms  can  be  used  in  an  OFDM  system  design  to 
induce  diversity  that  can  be  exploited  at  the  receiver  end. 

To  overcome  the  problems  of  a  simple  time  delay  diversity  scheme  in  a  frequency 
selective  fading  channel  for  MIMO-OFDM  systems,  a  new  approach  of  cyclic  delay 
diversity  was  suggested  in  [6].  This  simple  scheme  does  not  require  any  additional  guard 
interval  with  an  increasing  number  of  transmit  antennas.  The  combination  of  cyclic  delay 
diversity  and  OFDM  in  MIMO  systems  has  been  referred  to  as  multicarrier  delay 
diversity  modulation  (MDDM)  in  this  thesis  and  this  can  be  considered  a  special  type  of 
space-time  coding.  For  a  frequency  selective  channel,  a  cyclic  guard  interval  (cyclic 
prefix)  of  duration  G  is  added  at  the  beginning  of  each  OFDM  symbol.  This  guard 
interval  is  greater  than  or  equal  to  the  maximum  channel  delay  M  to  mitigate  the 
intersymbol  interference  (ISI).  The  orthogonality  of  subcarriers  is  paramount  for  OFDM. 
The  cyclic  prefix  converts  a  linear  convolution  channel  into  a  circular  convolution 
channel  and  the  interference  from  the  previous  symbols  will  only  affect  the  guard 
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interval.  This  restores  the  orthogonality  at  the  receiver.  Adding  zeros  as  the  guard  interval 
can  alleviate  the  interference  between  OFDM  symbols.  For  a  flat  fading  channel,  the 
guard  interval  can  be  eliminated  to  increase  the  data  rate.  [2,  22] 

The  space  time  code  for  MDDM  can  be  represented  as  [2] 

/  rp  rT  rF 

J'N-G  XN-G+ 1  •”  XN-1 

rp  rp  rp 

J'N~G- 1  ^ N-G  •”  -°N-2 

XN-G-L+ 1  XN-G+ 2  •••  XN-L 

v  cyclic  guard  interval 

where  G  is  the  length  of  the  guard  interval  and  N  is  the  number  of  points  in  the  IDFT. 

The  MDDM  scheme  has  some  advantages  over  space  time  block  coding  (STBC) 
and  space  time  trellis  coding  (STTC).  STBCs  require  all  block  codes  be  orthogonal  and 
full  rate  transmission  is  not  achievable  for  more  than  two  transmit  antennas.  The 
complexity  of  STTC  increases  exponentially  with  the  number  of  transmit  antennas. 
STBCs  and  STTCs  are  fixed  for  a  combination  of  transmit  and  receive  antennas  and  they 
do  not  offer  any  flexibility  in  changing  the  number  of  antennas  at  either  end.  For  each 
combination  of  the  number  of  antennas,  a  new  space-time  code  is  required.  These  codes 
were  designed  for  flat  fading  channels.  For  frequency  selective  channels,  they  cannot  be 
used  without  equalization  and  other  compensation  techniques.  However,  the  MDDM 
space  time  code  can  be  used  with  any  combination  of  transmit  and  receive  antennas  with 
little  modification  in  the  receiver  design.  The  MDDM  is  based  upon  the  OFDM 
modulation  scheme.  Therefore,  orthogonal  frequency  division  multiplexing  fundamentals 
are  discussed  very  briefly  in  the  following  sections.  [1,2] 

D.  ORTHOGONAL  FREQUENCY  DIVISION  MULTIPLEXING 

In  conventional  sequential  data  transmission  over  a  single  carrier,  the  data  rate  has 
a  limit  imposed  by  the  availability  of  frequency  spectrum  and  the  delay  spread  of  the 
channel.  The  bandwidth  required  is  inversely  proportional  to  the  data  symbol  duration, 
which  means  that  the  highest  achievable  symbol  rate  is  limited  by  the  available 

bandwidth  [16].  Considering  a  fixed  delay  spread  of  the  channel,  the  ISI  also  increases 

15 


Xn 


X 


N- 1 


X, 


Xn 


Xn 


X , 


XN~L+ 1  XN-L+2  XN-L+3 


data  block 


XN-1  ^ 


X 


N-2 


X 


N-L 


(2.32) 
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with  the  increase  in  the  symbol  rate  as  delayed  copies  of  the  symbols  coming  from  the 
multipath  can  have  significant  overlap  with  the  original  symbol.  The  lessening  of  ISI  will 
require  equalization  which  further  adds  to  the  complexity  of  the  system  [17].  These 
problems  can  be  mitigated  by  transmitting  the  data  in  parallel  on  multiple  carriers  with  a 
reduced  data  rate  on  each  carrier  compared  to  the  overall  data  rate.  The  reduced  data  rate 
will  require  reduced  bandwidth  which  should  not  be  more  than  the  coherence  bandwidth 
of  the  channel  to  avoid  frequency  selective  fading.  This  multicarrier  system  can  be 
designed  by  classical  frequency  division  multiplexing  (FDM)  [1].  In  this  scheme,  the 
carriers  need  to  be  well  apart  in  frequency  domain  to  avoid  inter-carrier  interference. 
Therefore,  a  guard  frequency  band  is  required  between  two  consecutive  subcarriers  which 
makes  this  scheme  highly  inefficient  in  frequency  spectrum  utilization.  This  problem  can 
be  eliminated  by  using  minimum-spaced  orthogonal  carriers.  In  OFDM,  all  carriers  are 
allowed  to  overlap  by  maintaining  orthogonality  of  all  the  subcarriers,  which  increases 
bandwidth  efficiency  [14].  The  comparison  of  bandwidth  efficiency  for  FDM  and  OFDM 
is  illustrated  in  Figure  5. 


bandwidth 


Ch.  1  Ch.  5 

/TYWYA 

◄ - 

bandwidth  savings 


(a) 


(b) 


Figure  5  Frequency  Spectrum  of  (a)  FDM  vs  (b)  OFDM  (After  Ref.  [14]) 

In  OFDM,  to  maintain  the  orthogonality  of  the  subcarrier  channels,  the  correlation 
between  signals  transmitted  on  subcarriers  must  be  zero.  Assume  that  the  available 
bandwidth  for  the  OFDM  system  is  A  W  and  it  is  to  be  divided  in  K  subcarriers.  The 
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input  serial  data  stream  is  to  be  converted  into  K  parallel  data  streams  which  are 
assigned  to  the  K  subcarriers  [1].  The  symbol  duration  of  the  input  serial  data  is  T  with 

serial  data  rate  of  f  s  =  VL<  .  Therefore,  if  the  number  of  parallel  data  streams  is  equal  to 

/  s 

the  number  of  OFDM  subcarriers,  the  symbol  duration  for  OFDM  will  be 

Ts  =  KTs.  (2.33) 

Equation  (2.33)  indicates  that  the  symbol  duration  of  an  OFDM  signal  is  K  times  larger 
than  that  of  single  serial  stream  symbol  duration.  Therefore,  the  OFDM  scheme  has  the 
inherited  advantage  over  single  carrier  modulation  techniques  to  mitigate  ISI  and 
frequency  selectivity  of  the  channel.  The  OFDM  transmitted  signal  S(t )  can  be  written 
as 


K-l 


S(t)  =  Y,A{dim  cosi27tfJ )  -  dQm  sin(2 nfmt)} 


m= 0 


K-l 

m=0 


dT 


dl  +  dQm 


^cos(2 nfj)  - 


dl  +  dQm 


^sin(2  nfj) 


(2.34) 


where  A  is  a  constant,  dj  ,  dQ  are  the  information-bearing  components  of  the  signal 
and  am  =  A^JVd^J) . 

Using  the  trigonometric  identity 

cos(a  +  P)  =  cos(a)  cos  (ft)  -  sin(«)  sin(/?)  (2.35) 

S(t)  can  be  written  as 


Sit)  =  X°mcos(2  7lfmt  +  0m) 

m= 0 


=  Yjam{COSi2^Lt)COSiGm) 

m= 0 


sm(2?F  fmt)  sin(6*m)} 


(2.36) 
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where  Q„.  =  tan  1 


.  The  correlation  between  any  two  symbols  transmitted  on 


separate  subcarriers,  represented  as  Rt] ,  must  be  equal  to  zero  to  maintain  the 
orthogonality  of  subcarriers  [1]. 

/•  +oo 

RtJ  =  L  s^s^dt 

=  £  ai  cos(2 +  0,t)ai  cos(2;r/;i  +  03t)dt 

Using  the  trigonometric  identity 


(2.37) 


cos(a)cos(/?)  =  ^-cos(a  +  /?)  +  ^-cos(«  -  /?) . 


(2.38) 


Equation  (2.37)  can  be  written  as 

-  ]£  (cos(2 +  f3)t  +  {6t  +  03)t)  +  cos(2 n{fi  -  f^t  +  {6t  -  d:i)t))lt  (2.39) 


R  _  a«°j  rT. 

ij  2  Jo 


where  an  a  ,  q,  and  q  ,  are  constant  for  the  symbol  duration. 


For  2 7i(ft  +  /,)»—,  Equation  (2.39)  can  be  written  as 


R:,  =  cos(2 -  fyt  +  [et  -  d:i)tyu . 


(2.40) 


From  Equation  (2.40),  Rt1  =  0  if 


(f,  ~  fjWs  =  M  where  M  e  set  of  positive  integers 
=•  f,-f=M/T 


(2.41) 


Therefore,  minimum  frequency  separation  between  two  consecutive  subcarriers  to 
maintain  orthogonality  must  be 


(2.42) 


where  fs  is  the  rate  of  OFDM  symbols. 
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1.  Generation  of  OFDM 

From  the  earlier  discussion  of  orthogonality  in  an  OFDM  signal,  all  subcarriers 
are  orthogonal  to  one  another  and  the  center  frequencies  of  any  two  subcarriers  differ  by 
a  positive  integer  multiple  of  a/.  In  spite  of  the  overlapping  of  subcarrier  frequency 
bands,  the  transmitted  symbol  can  be  recovered  at  the  receiver  end  without  any 
interference  due  to  zero  correlation  of  these  subcarriers.  In  OFDM,  serial  data  symbols 
are  converted  to  parallel  to  form  one  OFDM  symbol  to  be  assigned  to  K  subcarriers. 
Now,  the  bandpass  signal  can  be  represented  as 

S(t)  =  TJA  {d,m  cos (wj)  -  dQm  sin(wmf)}  (2.43) 

m= 0 

where  the  subcarrier  frequencies  wm  =  2 nfm  and  fm  =  f0  +  mi\f .  The  equation  (2.43) 
can  be  written  as 


K- 1 


iw^t 


S(t)  =  R  e\YjA{dIa  +idQm) 

tm=0 

=  R  eJ^dimle^l 

lm=0  J 

r  k-i 

=  Re 

lm= 0 
f  K- 1 

=  Re  ^d[m]ei2?r,iiV2^ 


^171= 0 

rI<-l 


=  R  e^(t)( 


i2nf0t 


<m=0 


(2.44) 


where  d[m]  =  AdI  +  iAdQ  and  the  complex  envelope  of  the  transmitted  signal  is 
denoted  by 

D(t)  =  Yjd[mV27[mAft.  (2.45) 

m= 0 

Matched  filters  or  correlation  demodulators  can  be  used  at  the  receiver  end  to 
recover  the  symbol  for  each  subcarrier.  Subsequently,  for  the  implementation  of  this 
scheme,  K  modulators  and  K  matched  filters  are  required.  As  the  number  of  subcarriers 
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increases,  the  complexity  of  the  system  also  increases.  For  a  large  number  of  subcarriers, 
the  complexity  of  the  system  makes  it  impractical.  [1] 


For  making  this  scheme  more  practical,  consider  sampling  the  complex  envelope 
signal  D  at  the  rate  of  K  /  Ts,  i.e.,  the  sampling  interval  is  Tj  K  .  Now  using  Equation 
(2.42),  Equation  (2.45)  can  be  represented  as 


D 


MTS 

K 


A'-l 


YJd\m]e‘2mnAfMrJK 


m= 0 


=  ^d[m]es"T‘/7yf 

m= 0 


K- 1 


m= 0 


(2.46) 


Comparison  of  Equations  (2.30)  and  (2.46)  clearly  yields  that  the  right  side  of  Equation 
(2.46)  is  the  IDFT  of  d[m\  scaled  by  K  . 


D  =  K IDFT  [d[m]}  (2.47) 

Equation  (2.47)  facilitates  the  implementation  of  an  OFDM  system  by  taking  the 
IDFT  of  the  original  data  stream.  Therefore,  implementation  of  the  OFDM  scheme  is 
completely  digital.  The  complexity  of  the  system  is  decreased  as  compared  to  the  multi¬ 
oscillator  based  OFDM  modulation  technique.  Even  for  large  76,  IDFT/DFT  can  be 
efficiently  implemented  by  using  fast  Fourier  transforms  (FFT)  which  make  computation 
much  faster.  [1] 

Similarly,  the  receiver  performs  the  FFT  operation  on  the  received  signal  to 
recover  the  original  parallel  data.  Employing  the  FFT  in  the  receiver  drastically  reduces 
the  complexity  as  compared  to  employing  coherent  demodulators.  The  receiver  can 
separate  the  subcarriers  simply  by  using  baseband  signal  processing  techniques.  These 
subcarrier  signals  can  be  integrated  for  the  duration  of  the  symbol  period  to  produce 
decision  variables  for  the  estimation  of  the  transmitted  data. 

2.  Cyclic  Guard  Interval 

As  previously  discussed,  the  OFDM  modulation  technique  has  the  advantage  of 
reducing  the  ISI  caused  by  a  multipath  channel.  This  is  achieved  by  converting  serial  data 
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at  a  high  rate  into  parallel  streams  each  at  a  lower  data  rate  and  increasing  the  symbol 
transmission  duration  on  orthogonal  carriers.  In  order  to  eliminate  the  ISI  completely  due 
to  the  time  delay  spread  of  the  multipath  channel,  a  guard  interval  is  added  before  each 
OFDM  symbol.  The  length  of  the  guard  interval  is  to  be  greater  than  the  anticipated  delay 
spread  of  the  channel  [1].  This  guard  interval  insertion  costs  both  in  terms  of  power  and 
data  rate  or  bandwidth.  If  the  length  of  the  guard  interval  is  less  than  the  delay  spread  of 
the  channel,  then  multiple  delayed  versions  will  induce  ISI  and  inter-carrier  interference 
by  affecting  the  orthogonality  of  the  subcarriers  due  to  overlapping  of  different  OFDM 
symbols  [14].  Taking  advantage  of  the  cyclic  nature  of  the  DFT  represented  in  Equation 

(2.31) ,  the  guard  interval  can  be  constructed  by  some  cyclically  shifted  portion  of  the 
OFDM  symbol.  Therefore,  it  is  referred  to  as  a  prefix  in  the  literature.  The  addition  of 
this  cyclic  guard  interval  has  already  been  depicted  in  the  space  time  code  (Equation 

(2.32) ).  Additionally,  this  cyclic  prefix  also  facilitates  the  synchronization  of  the  carrier 
frequency  and  timing  to  obtain  the  orthogonality  of  the  subcarriers  [12].  The  addition  of 
the  cyclic  guard  interval  is  illustrated  in  Figure  6. 
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Figure  6  OFDM  Symbols  with  Guard  Intervals 

E.  THE  MULTIPATH  AND  FADING  CHANNEL 

In  terrestrial  wireless  communications,  the  line  of  sight  (LOS)  path  for  signal 
reception  is  often  not  available.  The  signal  arriving  at  the  receive  antenna  is  comprised 
of  many  replicas  of  the  signal  coming  from  a  number  of  reflectors  and  scatters  present  in 
the  medium  of  transmission.  This  phenomenon  is  known  as  multipath.  Even  when  the 
LOS  path  signal  is  present,  the  multipath  components  are  usually  also  present.  [15,  16] 

In  this  section,  a  general  view  of  the  multipath  channel  is  presented  to  understand 
its  nature  and  a  baseband  equivalent  model  is  derived  for  use  in  simulation.  The  MDDM 
technique  will  be  tested  in  a  multipath  faded  channel  for  performance  analysis  and 
evaluation. 
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In  a  multipath  channel,  the  signal  travels  via  several  different  paths  before 
arriving  at  the  receive  antenna.  Due  to  the  different  path  lengths  involved,  the  amplitude 
and  phase  of  these  received  replicas  are  not  the  same.  For  example,  if  a  short  duration 
pulse  is  transmitted  in  a  multipath  channel,  then  a  train  of  pulses  of  different  amplitudes, 
phases  and  arrival  times,  will  be  received.  As  a  result,  the  received  signal  can 
significantly  vary  in  amplitude  and  phase  and  the  spectral  components  of  the  signal  are 
affected  differently  by  the  multipath  faded  channel.  Therefore,  the  frequency  response  of 
the  channel  may  not  be  flat  over  the  entire  bandwidth  of  the  signal.  Multipath 
components  of  the  signal  create  time-spreading  of  the  signal  and  cause  intersymbol 
interference.  [15,  16] 

The  channel  conditions  may  not  remain  constant  over  time  and  the  media 
composition  (e.g.,  stratosphere,  ionosphere)  may  also  change.  Furthennore,  the  number 
and  position  of  the  reflectors  and  scatters  cannot  be  considered  fixed.  In  mobile 
communications,  due  to  the  motion  of  transmitter  and  receiver,  the  multipath 
arrangements  can  never  be  assumed  constant  and  Doppler  shift  in  signal  frequency 
proportional  to  the  relative  velocity  is  also  observed.  Therefore,  the  multipath  fading 
channel  is  a  time  varying  channel  [15].  Time  varying  and  time  spreading  aspects  of  the 
channel  cannot  be  predicted  or  calculated  precisely.  For  creating  a  good  model  of  the 
channel,  these  parameters  are  measured,  and  based  on  these  measured  statistics,  the 
channel  can  be  characterized. 

Doppler  spread  in  the  signal  frequency  describes  the  time  varying  nature  of  the 
channel.  The  coherence  time  is  inversely  proportional  to  the  Doppler  shift.  The  coherence 
time  is  the  statistical  measure  over  which  the  channel  response  does  not  change 
significantly.  To  characterize  the  channel,  the  coherence  time  is  compared  with  the 
symbol  duration.  If  the  coherence  time  is  less  than  the  symbol  durations,  then  the  channel 
is  called  a  fast  fading  channel.  In  the  case  of  coherence  time  greater  than  the  symbol 
duration,  it  is  described  as  a  slow  fading  channel.  [15,  16] 

The  multipath  delay  spread  is  characterized  by  the  coherence  bandwidth. 
Coherence  bandwidth  has  been  defined  in  [17]  as,  “a  statistical  measure  over  which  the 
frequency  response  of  the  channel  is  considered  flat.”  If  the  bandwidth  of  the  signal  is 
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greater  than  the  coherence  bandwidth  of  the  channel,  then  different  frequency 
components  will  be  attenuated  differently  and  the  phase  variations  will  also  be  nonlinear. 
This  type  of  channel  is  said  to  be  a  frequency  selective  fading  channel.  If  the  coherence 
bandwidth  of  the  channel  is  greater  than  the  signal  bandwidth,  then  all  the  frequency 
components  will  face  flat  fading  with  almost  linear  phase  changes.  This  type  of  channel 
is  said  to  be  a  frequency  nonselective  fading  channel.  [15,  16] 


This  fading  channel  can  be  modeled  statistically  as  there  are  a  large  number  of 
variables  affecting  the  channel  response.  Most  of  these  variables  are  random  in  nature 
and  several  probability  distributions  can  be  considered  to  model  these  random  variables 
[16].  When  there  are  a  large  number  of  reflectors  and  scatters  in  the  physical  channel,  the 
cumulative  effect  of  this  large  number  of  random  variables,  as  per  the  central  limit 
theorem,  leads  to  a  Gaussian  process  model  for  the  channel  response.  In  the  case  of  no 
LOS  component  present,  the  process  has  zero  mean  with  magnitude  following  the 
Rayleigh  probability  distribution  and  the  phase  is  uniformly  distributed  on  [-  n  to  n\ . 
[15,  16] 

1.  Flat  Rayleigh  Fading  Channel 

A  slow  fading  frequency  nonselective  channel  was  simulated  to  test  the 
multicarrier  delay  diversity  modulation  scheme.  The  channel  fading  gain  was  kept  fixed 
for  the  symbol  duration  to  make  it  a  slow  fading  channel.  Considering  the  no  line  of  sight 
(LOS)  path  case,  a  flat  fading  channel  is  usually  simulated  using  a  Rayleigh  distribution 
for  the  magnitude  of  the  channel  response  [16].  The  Rayleigh  distribution  is  a  special 
case  of  the  Ricean  distribution  with  no  line  of  sight  component.  The  probability  density 
function  for  a  Rayleigh  random  variable  can  be  derived  from  the  Ricean  probability 
density  function  as  defined  in  [15] 


ac 

''  pvn 

r-(«c2  +«2P 

T 

(  aac) 

2  eXL 

a 

l  2a2  J 

10 

U2  J 

u{ac) 


(2.48) 


where  J0(.)  is  the  modified  Bessel  function  of  the  first  kind  of  zero  order,  u{.)  is  the 

unit  step  function  and  a 1  is  the  power  in  the  LOS  signal  component  and  the  average 
received  signal  power  is 
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s\t)  =  ac2  =  a2  +  2a2 


(2.49) 


jl  ac>0 

jo  ac  <  0. 


(2.50) 


For  the  Rayleigh  distribution  case,  there  is  no  line  of  sight  component  so  a  =  0  and 
/0(0)  =  1.  Using  Equation  (2.48),  the  Rayleigh  probability  density  function  can  be 
represented  as 


/a(«c)  =  ““ |"eXP 


<J 


f  2  \ 

~«c 

v2a2y 


u{ac)  . 


(2.51) 


For  the  baseband  simulation  of  a  Rayleigh  random  variable,  two  zero  mean 
independent  real  Gaussian  random  variables  were  summed  as  X  +  jY  [17].  The 
magnitude  of  this  complex  random  quantity  is  the  desired  Rayleigh  random  variable  and 
simulates  the  magnitude  of  the  channel  frequency  response.  The  derivation  of  the  proof 
that  the  complex  sum  of  two  zero  mean  Gaussian  random  variables  has  magnitude  with 
Rayleigh  distribution  and  phase  unifonnly  distributed  in  [-zr,  n]  is  largely  based  on  [8], 


Let  X  and  Y  be  two  zero  mean  independent  identically  distributed  (IID)  Gaussian 
random  variables  and  their  complex  sum  is  represented  by 

Z  =  X  +  iY.  (2.52) 

As  Z  is  to  simulate  the  frequency  response  of  a  Rayleigh  flat  fading  channel,  it  can  be 
written  as 


Z  =  hew  (2.53) 

V  =  \zf  =  X'2  +  Y'2 .  (2.54) 

Since  X  and  Y  are  zero  mean  IID  Gaussian  random  variables  (GRVs) 

X  =  Y  =  0  (2.55) 

<j2x  =  a2  =  a2 .  (2.56) 
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The  random  variable  V  is  a  sum  of  two  squared  zero  mean  GRVs  and  is  called  a  central 
chi-squared  random  variable  of  degree  2.  The  probability  density  function  for  a  central 
chi-squared  random  variable  of  degree  n  can  be  given  as  [16] 


fv  (y) 


_ 11  w/2  T 

(J  z  1 


v(,/2)-Vv/2  <rM(v) 


(2.57) 


w 

where  n  is  the  number  of  independent  variables  and  T(.)  is  the  gamma  function.  For  a 

central  chi-squared  random  variable  of  degree  two  (n  =  2),  Equation  (2.57)  can  be 
rewritten  as 


/r  (v) 


2cr2r(l) 


e~vl2a~u(v). 


(2.58) 


Since 


r(i)  =  (i— i)!=i , 


(2.59) 


[15],  the  probability  density  function  of  V  as  per  Equation  (2.58)  can  be  written  as 

fv{v)  =  ^-re ~V  l2a' «(v).  (2.60) 

Za 

The  magnitude  of  the  simulated  Rayleigh  flat  faded  channel  frequency  response  can  be 
written  as 

h=\Z\  =  y/V  (2.61) 

The  probability  density  function  for  h  can  be  obtained  by  transforming  the  probability 
density  function  of  V  given  in  Equation  (2.58)  according  to  [18] 

(2.62) 

v=h2 


/„«= 


i 


dh  /  dv 


Mv ) 


where 


dh  _  1 

~dv  =  2(v)1/2  ‘ 


(2.63) 
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Substituting  Equations  (2.60)  and  (2.63)  into  Equation  (2.62)  yields 


r  (l\  2\l/2  1  -h2 /2a2 

fH(h)  =  2(h)  —^e  u 

1(7 

(7 

Comparison  of  Equations  (2.51)  and  (2.64)  clearly  indicates  that  the  magnitude  of  the 
complex  sum  of  two  zero  mean  IID  GRVs  follows  the  Rayleigh  distribution.  The  mean 
and  variance  of  h  as  given  in  [18]  are 


« 

(2.64) 


(2.65) 


(2.66) 


Now,  the  probability  distribution  function  for  the  phase,  0 ,  of  Z  is  to  be  derived. 
The  phase  0  can  be  represented  as  [8] 

-if  Y 

tan  — 

u 

/  Y 
tan  — 

u 

/  Y 
tan  — 

u 

Let  the  ratio  of  zero  mean  IID  GRVs  X  and  Y  be  defined  as  [8,  19] 

R  =  —  .  (2.68) 

X 

The  joint  probability  density  function  of  X  and  Y  can  be  written  as 

fxr  (x,  y)  =  fx(x)fY  (■ y )  =  fx(x)fx  (y)  (2.69) 

where 


—  n  if-oo<X<0  and  -oo  <  Y  <  0  (referred  to  as  case  A) 

if  0  <  X  <  oo  and  -oo  <  Y  <  oo  (referred  to  as  case  B)  (2.67) 

+  7T  if-oo<X<0  and  0  <  Y  <  oo  (referred  to  as  case  C) 


(2.70) 
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For  case  A,  -oo  <  X  <  0  and  -  oo  <  Y  <  0  which  implies  -n  <  9  <  -n  /  2 .  The 
conditional  cumulative  distribution  function  of  R  can  be  represented  as  [8,  19] 


Fr\a  {r )  =  Pr  {R  <  r  [ A}  =  Pr  <  H- A 

=  \0_Sr/XY\^X,y)dydX 

=  \r  fXY{x,y)dydx 


(2.71) 


where  it  is  noted  that 


fxY\A  (A’T) 


4  fxY  {x,y)  if  -  oo  <  x  <  0  and  -  oo  <  y  <  0 


otherwise 


(2.72) 


The  conditional  probability  density  function  can  be  obtained  by  taking  the  derivative  of 
the  conditional  cumulative  distribution  function  with  respect  to  random  variable  r  [8,  19] 


fR  iaW  = 


dFR\A  (r)  Ad  ffO 


//  /*0  fO  ,  . 

4—  f  f  fxr  (x,y)dydx  . 
dr  LJ“°oJ,'r  J 


From  Leibniz’s  rule  [8,  19] 


fR |A  ir)  =  4J_xi  -%fxY  rx) dx  ■ 


Substituting  Equation  (2.69)  and  (2.70)  into  Equation  (2.74)  yields  [8,  19] 


1  r° 


A(0=-— r  1 


(2x}dx  . 


(2.73) 


(2.74) 


(2.75) 


Equation  (2.75)  can  be  represented  as 


For  case  A, 


0  =  tan  1  ( r}-K  , 

d9_  1 

dr  1  +  r2 
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(2.76) 


(2.77) 


(2.78) 


The  conditional  probability  density  function  of  0  can  be  attained  from  Equation  (2.76) 


/o|A  (0)  =!,,„/  ,?„|/Ji|A(r) 


dO  /  dr 


Jr  |a  v 


=tan(0) 


4|.W  =  (l  +  ''2)-7— IT 
nt  l  +  r  1 

Equation  (2.80)  can  be  written  as 


(2.79) 


(2.80) 


-=tan(0) 


/«*(»)=- 

n 


(2.81) 


which  implies 


hw = A 

Zn 


for  -n  <  9  <  — . 


(2.82) 


Similarly  forO  <  X  <  oo,-co  <  7  <  oo  and  for-oo  <  X  <  0,-qo  <  Y  <  0 ,  it  can  be 
proven  [8]  that 


U(o)=- 

n 

(2.83) 

n 

(2.84) 

Therefore,  the  unconditional  probability  density  function  for  the  phase  6  can  be  written 
as 

/eM  =  /0|A(^)Pr(A)  +  /0|B(^)Pr(B)  +  /0|C(^)Pr(C) 
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0  otherwise 
where  /(.)  is  the  indicator  function  defined  as 
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1  if  x  e  A 

0  otherwise 


(2.86) 


In  this  section,  it  was  proved  that  the  complex  sum  of  two  zero  mean  IID  GRVs 
gives  a  complex  sum  with  Rayleigh  distributed  magnitude  and  unifonnly  distributed 
phase  in  (-zr,  n] .  Therefore,  this  model  can  be  used  to  simulate  the  frequency  response 

of  a  multipath  flat  Rayleigh  faded  channel. 

2,  Maximal-Ratio  Combining 

As  stated  in  previous  sections,  many  wireless  communication  systems  operate  in 
multipath  channels  and  the  performance  in  a  multipath  channel  is  often  reduced  as 
compared  to  an  AWGN  channel.  Diversity  techniques  can  be  employed  to  mitigate  the 
effect  of  multipath.  Diversity  simply  implies  the  transmission  or  reception  of  multiple 
copies  of  the  same  signal.  Diversity  can  be  achieved  in  time,  frequency  and  space 
domains.  It  is  assumed  that  all  these  diversity  receptions  are  independent  each  with  an 
independent  channel  response .  [15] 

For  BPSK,  the  transmitted  baseband  symbols  are  defined  as 


x0  (V)  =  Ae°  =  A  for  0  <  t  <  Tc 

xl  (t)  =  Ae~m  =-A  for  0  <  t  <  Tc 


(2.87) 


where  Tc  is  the  symbol  duration  in  each  time  diversity  and  A  is  the  peak  amplitude  of 
the  BPSK  signal.  The  energy  for  each  time  diversity  reception,  E ,  can  be  given  as  [16] 

E  =  J7"1  \xk(t)f  dt  =  A2Tc  for  k  =  0, 1 .  (2.88) 

For  time  diversity,  it  is  assumed  that  each  received  diversity  signal  passes  through 
an  independent  faded  channel  and  can  be  represented  as 

rl(t2)  =  hle'e'xk(t)  +  nl(t)  for  0<t  <Tc,  k  =  0,\  (2.89) 

where  l  represents  the  number  of  the  diversity  reception  and  n,(t)  is  the  complex  valued 
AWGN  with  the  circularly  symmetric  probability  density  function.  The  power  spectral 
density  function  for  nt(t)  is 

Snn(f)  =  Nt 
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(2.90) 


where  the  power  spectral  density  functions  for  the  real  and  imaginary  components  of 
n,  (t)  are 


^Re[n]Re[n]  (/)  ‘^Im[n]Im[n]  (/)  ' 

For  all  the  diversity  receptions,  the  received  signal  can  be  written  as 

r(t)  =  X[^e^%(f)  +  n;(t)]. 

/=! 


(2.91) 


(2.92) 


The  random  variable  Ylk  after  the  correlation  receiver  can  be  represented  as  [8] 


Yi,k  =A\0  ( h,e -,e'  xk  ( t )  +  n,  (t))dt 

jA2Tch,e~Wl  +N,  for  k  =  0  (2-93) 

\-A2Tch,e-Wl  +N,  for  k  =  1, 


where  Nl  is  a  zero-mean  complex  Gaussian  random  variable  with  circularly  symmetric 
probability  density  function  which  represents  the  noise  component  and  can  be  given  as 

Nl  =  n,  ( t )  dt.  (2.94) 

Since  the  integrator  is  a  fdter  with  frequency  response  Hint(f )  and  impulse  response 


fl  if  0  <  t<Tb 
1 0  otherwise. 


(2.95) 


Then,  the  variance  of  can  be  calculated  using  Paresval’s  Theorem  and  Equations 
(2.90)  and  (2.95) 

+oo 

<  =  e{|«,2|}  =  a2  1 1 hm  (/)|2  sm  (f)df 

-00 

+00 

-Jl  I  hlMN«df  ■  (2-%) 

-00 

=  Th  =  EN0 


Substituting  Equation  (2.88)  into  Equation  (2.93)  yields 
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for  k  =  0 
for  k  =  1 . 


(2.97) 


{Ehfi-*1  +N, 
\-Eh,eiei  +N, 


It  is  assumed  that  the  exact  channel  state  information  (CSI)  is  known  at  the 
receiver  end.  The  complex  conjugate  of  the  CSI  is  multiplied  by  the  random  variables  Y,  k 

to  produce  the  random  variables  Z,  k .  Thus,  the  phase  shift  in  the  channel  is  compensated 
and  the  value  of  the  mean  of  the  random  variables  Zlk  is  proportional  to  the  signal 

power.  Therefore,  a  strong  received  signal  carries  a  larger  weight  than  a  weak  received 
signal.  All  the  diversity  receptions  are  added  to  form  the  random  variable  Zk .  This 
optimum  combiner  is  called  the  maximal  ratio  combiner  (MRC)  [16].  Figure  7  illustrates 
the  MRC  for  BPSK  with  time  diversity. 


Combiner 


Figure  7  MRC  for  BPSK  with  Time  Diversity  (After  Ref.  [8,  15]) 


After  maximal  ratio  combining,  the  random  variable  Zk  can  be  written  as 


1=1 

-£[£/>=  +  A/>AT,] 


for  k  =  0 

for  k  =  1 . 


(2.98) 


The  real  part  of  Zk,  denoted  by  £k  =  Re  { Zk } ,  is  fed  to  the  comparator  for  comparison 
with  the  predefined  threshold  level.  For  BPSK,  if  £k  >  0,  the  receiver  decides  k  =  0  was 
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transmitted.  If  £k  <  0,  the  receiver  decides  k  =  1  was  transmitted.  The  decision  variable 
£k  can  be  written  as 


C=Re[Z,]  = 


L 

+Re[/*/e,'*'A/]) 

i=i v 

for  k  =  0 

( Ehf  +  Re  [  h,ew'  7V;  ]  ) 

for  k  =  1 

(2.99) 


1=1 


From  Equation  (2.99),  it  is  evident  that  the  decision  variable  £k  is  a  Gaussian  random 

L 

variable  with  conditional  mean  (conditioned  on  the  value  of  ^  hf  ) 

1=1 

S  =  (2-100) 


1=1 


Since  the  N,’ s  are  independent  zero-mean  complex  Gaussian  random  variables  with 
circularly  symmetric  probability  density  function  and  variance  equal  to  EN0 ,  the 
conditional  variance  of  £k  can  be  expressed  as: 


2  2  2 

°V  =  °V0  =  GCx 


=  ^  Var  |  Re  ( h^1  N, ) 

/=i 

=  Z*,2Var{Re(JV,)}  . 


^tyVarjAf,} 

^  /=  1 


ENn 


Yjh 


i=\ 


(2.101) 


Assuming  the  probability  of  transmitting  a  “1”  bit  and  a  “0”  bit  are  equal  and  using  the 
fact  that  the  conditional  probabilities  of  bit  error  Phk  are  equal  due  to  the  symmetry  of  the 

noise  probability  density  function  and  the  zero  threshold,  it  is  possible  to  calculate  the 
conditional  probability  of  bit  error  as: 
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n=^w+^ni,=^«=p>-{f.<o} 


=  Pr 


=  Pr 


=  0 


=  0 


u  ^  u 

<7, 

Vo" 

■Co 

>u 

fT, 

-c\ 

7  \ 

ac  ) 

- X 

(2.102) 


Ifl'v 

7Vo  /= 1 


Where  Q  (x)  is  defined  as  [18] 


Let  y  be  defined  as 


Q{x)=^y~2du 


E  Cx  ,2 

r=Y^h'  • 
iV0  1=1 


(2.103) 


(2.104) 


Substituting  Equation  (2.103)  into  (2.102),  it  is  possible  to  rewrite  the  conditional 
probability  of  bit  error  as  a  probability  of  bit  error  conditioned  on  y : 


P*(r)=Q(-j2r )• 


(2.105) 


The  average  probability  of  bit  error  can  be  obtained  by  taking  the  expectation  of  Pb  (y) 
with  respect  to  random  variable  y  [15,  16].  The  average  probability  of  bit  error  can  be 
written  as 


n=Epito] 

=  £P‘b)/r(j'Vr 


(2.106) 

where  fr(y)  is  the  probability  density  function  for  y .  Similarly,  yt  can  be  written  as 
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(2.107) 


The  probability  distribution  function  of  hf  was  derived  in  the  last  section.  Taking///  =  vt. 
Equation  (2.60)  can  be  represented  as 


(2.108) 


The  average  SNR  per  diversity  reception  can  be  written  as 


Eh 2 

r,=E[r,]=  e  ^ 


(2.109) 


The  expectation  can  be  evaluated  as 


E[vi]  =  \vifvSvi)dvi 


f  1 

=J  V,r^ 


-e-v'n^dvl=2cj2. 


(2.110) 


Now  the  Equation  (2.109)  can  be  written  as 


_  IEcj2 

ri=^r~ 


(2.111) 


The  characteristic  function  of  Vl  can  be  written  as  [8,  16] 


F,;W  (l-/2a  V)' 


(2.112) 


If  L  IID  random  variables  are  added,  then  the  probability  density  function  of  the  sum  is 
the  L  -fold  convolution  of  the  probability  density  function  of  the  single  random  variable. 
Therefore,  the  characteristic  function  of  the  sum  is  the  characteristic  function  of  the 
single  random  variable  raised  to  the  power  of  L  [16,  18].  Therefore,  the  characteristic 


L, 

function  for  V  =  ^  Vt  is 
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1 


(2.113) 


Fv  ( <a ) 


(l-z'2cr7&>) 


The  probability  density  function  of  V  is  the  inverse  Fourier  transform  of  the 
characteristic  function  and  can  be  written  as  [8] 

L- 1 


fv  (v) = ■ 


-v!2<7* 


2L  <j2l  (L-\)\ 

Consistent  with  Equation  (2. 104),  it  is  possible  to  write 

Ev 


u(v) . 


r  = 


Nn 


(2.114) 


(2.115) 


and 

dy  _  E 
dv  N0 

Now,  the  probability  distribution  function  for  y  can  be  given  as 


(2.116) 


fr(r)  = 


-v/2<7 


^2lcj2l{L-\)\ 

Nn 


f(r)  = 


L- 1 


r 


r 

E2<t2 


f  E2cj2^L 

V^vT  y 


(Z-1)! 


Substituting  Equation  (2.111)  into  Equation  (2.117)  yields 


fr{r)  =  — 


7 


(r) 


-e  y,u 


W- 


(2.117) 


(2.118) 


Now,  substituting  Equation  (2.105)  and  (2.118)  into  Equation  (2.106),  the  average 
probability  of  bit  error  can  be  represented  as  [8] 
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(2.119) 


Ph  = 


(r,)L(L-\)\ 


_r_ 

e  7'dy . 


The  solution  for  Equation  (2. 1 19)  has  been  given  in  [16]  as 


Pb  = 


(l  —  ll) 

L  L- 1, 

fz-l+A 

2 

Zl 

1=0 

l  /  , 

(1  +  u) 
2 


where  u  is  defined  as 


(2.120) 


(2.121) 


In  time  diversity,  the  total  bit  energy  received  is  proportional  to  the  number  of  diversity 
receptions  and  is  represented  as 

Eb  =  LE  (2.122) 

where  L  is  the  number  of  diversity  reception  and  E  is  the  energy  per  diversity  reception. 
Substituting  Equation  (2.122)  into  Equation  (2.109),  yl  can  be  written  as 


r,  = 


2 

LN0 


(2.123) 


Similarly,  in  the  case  of  space  diversity  of  order  J  as  illustrated  in  Figure  8,  each 
diversity  reception  is  processed  by  a  separate  correlator.  The  integrating  time  period  in 
each  correlator  is  Tb .  Therefore,  the  bit  energy  per  diversity  reception  can  be  represented 
as 


E  =  Eb 


(2.124) 


Now,  the  average  SNR  per  diversity  reception  y,  can  be  written  as 


r,  = 


2Eha2 

N0 


(2.125) 
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Combiner 


3  =  1 


Figure  8  MRC  for  BPSK  for  Space  Diversity  (After  Ref  [8,  15]) 


In  this  thesis,  the  performance  of  MDDM  in  MISO  and  MIMO  systems  is 
compared  with  the  SISO  system.  In  the  case  of  the  SISO  system,  there  is  no  diversity  and 

L  =  1 .  Thus,  the  probability  of  bit  error  Pb  as  per  Equation  (2. 120)  is  given  as  [15,  16] 


Pb  = 


1  —  u 


Substituting  Equation  (2.121)  into  Equation  (2.126)  yields 

r  i  — 


Pb  = 


r, 


K  \l  +  Yl 


(2.126) 


(2.127) 


The  probability  of  bit  error  for  a  BPSK  MISO  system  in  a  multipath  fading  channel  as 
given  in  Equation  (2. 127)  will  be  compared  with  the  simulated  and  theoretical  probability 
of  bit  error  for  MDDM  MISO  and  MIMO  systems. 

F.  SUMMARY 

In  this  chapter,  the  fundamentals  of  SISO,  SIMO  and  MISO  systems  were 
introduced  to  facilitate  better  understanding  of  MIMO  systems.  To  approach  the 
maximum  capacity  of  MIMO  channels,  space  time  codes  can  be  designed  in  conjunction 
with  the  modulation  techniques  to  incorporate  diversity.  Orthogonal  frequency  division 
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multiplexing  (OFDM)  can  be  used  to  mitigate  the  frequency  selectivity  of  the  channel. 
OFDM  utilizes  the  discrete  Fourier  transfonn  (DFT)  of  the  baseband  signal.  The 
multicarrier  cyclic  delay  diversity  technique  can  be  implemented  by  exploiting  the  cyclic 
shift  property  of  the  DFT.  A  baseband  model  for  the  flat  Rayleigh  fading  channel  was 
derived  for  simulation.  Lastly,  the  MRC  as  an  optimum  diversity  combining  technique 
was  discussed.  The  next  chapter  discusses  the  design  of  the  simulated  MDDM 
transmitters  and  receivers. 
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III.  MULTICARRIER  DELAY  DIVERSITY  MODULATION 
TRANSMITTER  AND  RECEIVER  MODELS 


This  chapter  introduces  the  implementation  of  the  multicarrier  delay  diversity 
modulation  scheme  at  the  block  level.  The  basic  design  of  the  transmitter  and  receiver  is 
discussed  briefly.  In  this  thesis,  the  MDDM  scheme  is  implemented  using  BPSK 
modulation.  The  inclusion  of  forward  error  correction  (FEC)  and  higher  bandwidth 
efficient  modulation  schemes  can  be  considered  for  future  work.  The  next  chapter 
illustrates  the  simulation  and  analysis  of  the  MDDM  transmitter  and  receiver. 

A.  THE  MULTICARRIER  DELAY  DIVERSITY  MODULATION  SCHEME 

The  scheme  of  implementing  delay  diversity  with  OFDM  for  MIMO  systems  was 
recommended  for  flat  fading  channels  in  [5].  The  cyclic  delay  diversity  approach  with 
OFDM  was  suggested  for  the  frequency  selective  fading  channel  in  [6],  In  this  research, 
cyclic  delay  diversity  with  OFDM  is  investigated.  The  implementation  of  multicarrier 
delay  diversity  in  MIMO  systems  has  been  discussed  in  the  previous  chapter.  For 
modeling  purposes,  the  MIMO  configuration  of  two  transmit  and  two  receive  antennas  is 
discussed  in  this  chapter.  For  simulation  and  analysis,  the  number  of  receive  antennas 
varies  from  one  to  three.  To  facilitate  better  understanding  and  precise  representation  of 
signals  from  the  transmit  antenna  to  the  receive  antenna,  the  following  notations  will  be 
used: 

•  Xk  represents  the  BPSK  symbol  at  k  -th  interval  before  the  multicarrier 
delay  diversity  modulation  and  after  demodulation 

•  Xlk  represents  the  BPSK  symbol  at  k  -th  interval  from  the  transmitting 
antenna  /  after  multicarrier  delay  diversity  delay  modulation 

•  x'm  represents  the  signal  x  to  be  transmitted  from  transmit  antenna  l  in 
the  rnxil  transmission  interval, 

•  rjt  represents  the  received  signal  at  receive  antenna  j  at  the  m~th  time 
interval, 

•  hlj  represents  the  channel  response  from  transmit  antenna  l  to  receive 
antenna  j, 

•  nJ  represents  the  AWGN  at  receive  antenna  j ,  and 

•  asterisk  (  )*  represents  the  complex  conjugate. 
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B.  MDDM  TRANSMITTER 

The  multicarrier  delay  diversity  scheme  will  be  simulated  with  binary  phase  shift 
keying  modulation.  A  block  diagram  of  the  MDDM  transmitter  with  two  transmitting 
antennas  is  shown  in  Figure  9.  The  function  of  each  block  is  explained  briefly  in  the 
following  discussion. 


Tx.  Antenna  1 


Figure  9  Block  Diagram  of  MDDM  Transmitter  (After  Ref.  [2]) 

1.  Binary  Information  Source  and  M  PSK  Modulator 

The  binary  information  source  generates  bits  0  and  1  with  equal  probabilities 

Pr[0]  =  Pr[l]  =  1/2 .  (3.1) 

Then,  the  bits  are  mapped  to  I  and  Q  symbol  coordinates  by  the  M  PSK  modulator. 
The  equivalent  lowpass  (i.e.  complex  envelope)  M  PSK  symbols  can  be  represented  as 

Xk(t)  =  AejMm~1)/MPT(t)  for  m  =  1, 2, . . . , M  (3.2) 

where  A  is  the  amplitude  of  the  signal,  M  is  the  number  of  possible  phases  (for  BPSK 
M  =  2,  for  QPSK  M  =  4),  T  is  the  symbol  duration,  and  PT (t)  denotes 

fl  for  0  <  t  <  T 

=  {o  otherwise  '  (33) 

2.  OFDM  Modulator 

Mapped  symbols  are  fed  to  the  OFDM  modulator  block.  In  this  block,  first  the 
serial  input  data  stream  is  converted  to  parallel.  Then,  the  IFFT  operation  is  perfonned  to 
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realize  OFDM  as  discussed  in  the  last  chapter.  The  FFT  size  denotes  the  number  of 
subcarriers.  Figure  10  illustrates  the  OFDM  modulator.  [2] 


Figure  10  Orthogonal  Frequency  Division  Multiplexing  Modulator 


In  this  thesis,  the  number  of  subcarriers  is  256  to  simulate  OFDM,  in  accordance 
with  the  IEEE  802.16a  standard.  In  this  standard,  three  types  of  subcarriers  are  defined: 
data  carriers  for  information  symbol  transmission,  pilot  carriers  for  the  estimation  of  CSI 
and  guard  carriers  (null  carriers)  placed  on  both  sides  of  the  frequency  spectrum  to  avoid 
intercarrier  interference  from  neighboring  frequency  bands.  The  assignment  of  the 
subcarriers  in  this  thesis  has  been  adopted  from  the  IEEE  802.16a  standard  and  given  in 
Table  1.  Figure  1 1  illustrates  the  subcarriers  organization  at  the  input  of  the  IFFT  block. 


Size  of  FFT 

256 

#  of  information  subcarriers 

192 

#  of  pilot  subcarriers 

8 

#  of  null  subcarriers  (including  the  DC 
subcarrier) 

56 

#  of  lower  frequency  guard  subcarriers 

28 

#  of  higher  frequency  guard  subcarriers 

27 

Frequency  indices  of  null  subcarriers 
(including  the  DC  subcarrier) 

-128,-127,. ..,-101,0 

+101, +102,. ..,+127 

Frequency  indices  of  pilot  subcarriers 

1,2,..., 8 

Table  1  Assignment  of  OFDM  Subcarriers  (After  IEEE  802.16a  standard,  Ref.  [20]) 
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28  null  subcarriers 


1 92  information 
symbols 

DC  subcarrier. 


8  pilot  subcarriers 


27  null  subcarriers 


- 128 

-101 

-100 

-60 

-36 

-1 

0 

IFFT 

+9 

+36 

+60 

+  100 

+  101 

+  127 

Figure  1 1  Assignment  of  Subcarriers  at  the  Input  of  IFFT  Block  (After  Ref.  [20]) 

The  output  of  the  IFFT  block  can  be  represented  according  to  the  definition  of  Equation 
(2.30)  as 

i  JV-l 

xm  =-YJXke»*km/N  for  m  =  0,1, 2, ...,255  (3.4) 

A  k= 0 

where  N  =  256,  and  the  index  of  the  subcarrier  is  represented  from  0  to  255  instead  of 
-128  to  128  for  convenience.  After  the  IFFT  operation,  the  OFDM  sequence  is 
converted  from  parallel  to  serial.  [23] 

3.  Cyclic  Delay  Addition 

After  the  IFFT  operation,  the  symbols  xn  are  cyclically  shifted  to  realize  a  cyclic 
delay  of  one  information  symbol  period.  For  transmitting  antenna  1 ,  there  is  no  delay  and 
D  =  0.  For  each  subsequent  antenna,  the  cyclic  delay  is  increased  one  information 
symbol  period.  The  signal  after  cyclic  delay  is  given  as 

xln  =  X(n-D)N  where  l  =  1,2 and  D  =  1-1  (3.5) 

4.  Guard  Interval  (Cyclic  Prefix)  Addition 

As  discussed  in  the  previous  chapter,  a  guard  interval  at  the  beginning  of  each 
OFDM  symbol  is  added  to  mitigate  the  effect  of  the  multipath  channel.  The  length  of  the 


42 


guard  interval  TG,  in  the  fonn  of  a  cyclic  prefix,  is  chosen  to  be  greater  than  or  equal  to 
the  anticipated  delay  spread  of  the  channel.  The  addition  of  a  guard  interval  is  an 
overhead  at  the  cost  of  data  rate.  The  guard  length  is  specified  with  respect  to  the  OFDM 
symbol  period  Ts  and  is  also  referred  to  as  guard  ratio  T(;  j  T  .  The  addition  of  the  guard 
interval  has  been  illustrated  in  Figure  6.  For  example,  in  the  IEEE  802.16  wireless 
standard,  of  the  guard  ratios  is  chosen  to  be  1  /  32 .  The  FFT  block  size  is  256.  Therefore, 
the  guard  length  is  to  be  8.  The  cyclic  prefix  constitutes  the  last  eight  samples  of  the  IFFT 
output  and  it  is  concatenated  to  the  beginning  of  the  OFDM  symbol  sequence.  Now,  the 
sequence  with  added  guard  interval  for  each  transmitting  antenna  can  be  represented  as 

xLg)  =  4248]  ^[249]  ...x[255]  x[0]  x[l]  x[2]...x[255}} 

(3.6) 

Xl+G)  =  4247]  ^t248]  -4254]  x[255]  x[0]  x[l]  -4254]]. 

5.  Digital  to  Analog  Conversion  and  RF  Modulation 

The  IFFT  operation,  the  cyclic  delay  operation  and  the  guard  interval  operations 
are  performed  in  the  discrete  time  domain.  For  transmission  purposes,  these  discrete 
symbols  are  converted  to  an  analog  signal  using  a  digital  to  analog  (D/A)  converter.  After 
digital  to  analog  conversion,  these  symbols  are  continuous  time  baseband  OFDM 
symbols.  Then,  these  symbols  are  upconverted  to  the  transmission  radio  frequency.  After 
requisite  amplification  to  the  desired  power  level,  these  symbols  are  transmitted  from  the 
respective  antennas.  In  this  thesis,  the  simulation  is  performed  at  the  discrete  time 
baseband  level  and  no  digital  to  analog  or  RF  modulation  is  implemented. 

C.  MDDM  RECEIVER 

The  multicarrier  delay  diversity  scheme  can  be  employed  with  any  number  of 
receive  antenna  without  changing  the  transmitted  space  time  code.  In  this  section,  the 
receiver  design  for  two  receive  antennas  is  discussed.  It  is  assumed  that  the  signals 
received  at  both  receiving  antennas  are  independent  of  each  other  due  to  independent 
channel  responses  and  no  antenna  correlation.  The  MDDM  receiver  has  almost  the  same 
structure  as  the  MDDM  transmitter  but  the  operations  are  performed  in  reverse  order.  The 
MDDM  receiver  design  is  illustrated  in  Figure  12.  At  the  receiving  antenna,  the 
cyclically  time  delayed  OFDM  symbols  are  received  after  passing  through  the  channel 
with  the  AWGN  and  can  be  represented  as 
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for  n  =  0,1,2, ...,N  +  G  - 1 
for  n  =  0, 1,2,..., TV  +  G  -  1 


(3.7) 


r}i+r)  =  A(hn  *  x1  +  h 21  *  x2  )  +  n1 

ln+tO  \  (»+G)  (n+C)  / 

r2  =  A(h'2*x'  +h^*x2  )  +  n2 

(n+G)  ^  (b+g)  („+g)  J 


where  *  represents  circular  convolution  (due  to  the  circular  guard  interval).  [2] 

After  RF  demodulation  and  analog  to  digital  (A/D)  conversion,  the  OFDM 
symbol  received  by  each  receiving  antenna  is  in  the  discrete  time  domain.  The  guard 
interval  added  at  the  transmitter  is  removed.  Then,  this  signal  is  fed  to  the  OFDM 
demodulator  as  a  baseband  signal.  In  the  OFDM  demodulator  block,  the  received  OFDM 
symbol  is  converted  from  serial  to  parallel  for  the  FFT  operation.  The  OFDM 
demodulator  is  illustrated  in  Figure  13.  After  the  FFT  operation  as  defined  in  Equation 
(2.29),  the  signal  can  be  represented  as 

255  2 

R[  =  FFT {r;(}  =  J](2( hlJ  *  ®')  +  nJ)e^"  ™  for  n  =  0,1, 2, ...,255.  (3.8) 

n= 0  1=1 


Rx.  Antenna  1 


Figure  12  Block  Diagram  of  MDDM  Receiver 
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Figure  13  Orthogonal  Frequency  Division  Multiplexing  Demodulator 


Using  the  linearity  property  of  the  discrete  Fourier  transform,  Equation  (3.8)  can  be 
represented  as 

Rl  =  HnX\  +  H2lX\  +  N\  for  k  =  0, 1, 2, . . . ,  255 

(3.9) 

R2k  =  Hl2X\  +  H22X2k  +  N2k  for  k  =  0, 1, 2, . . . ,  255 


where  Hlj  is  the  frequency  response  of  the  channel  from  transmit  antenna  l  to  receive 
antenna  j  and  Nj  is  the  AWGN  at  receive  antenna  j  in  the  frequency  domain.  Using 
the  circular  time  shift  property  of  the  discrete  Fourier  transform  as  defined  in  Equation 
(2.31),  Equation  (3.9)  can  be  written  as 


Rl  =  HnX\  +  H21Xle~i2M/N  +  N1 
Rl  =  HX2X{  +  E22X\e-i2nkD'N  +  N2 
Rl  =  ( Hn  +  E21eankDIN)X\  +  N1 
Rl  =  ( H 12  +  E22e-a*kD/N)Xl  +  N2 


for  k  =  0,1, 2, ...,255 
for  k  =  0,1, 2, ...,255 
for  k  =  0,1, 2, ...,255 
for  k  =  0,1, 2,... ,25 


(3.10) 


where  D  =  1  is  the  cyclic  time  delay  induced  at  transmit  antenna  two.  Now,  the  effective 
channel  response  at  the  receive  antenna  j  is  defined  as 

H3k  =  Hi1  +  H2kje-i2*k/N  .  (3.11) 

Substituting  Equation  (3.11)  into  Equation  (3.10),  the  signal  at  each  receiving  antenna 
can  be  written  as 
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(3.12) 


R\  H\X\  +  N\ 

Rl=HlXl+Nl 

Equation  (3.12)  shows  clearly  the  space  diversity  receptions  at  the  receiving  end. 
They  are  combined  by  using  the  optimum  maximal  ratio  combining  receiver.  MRC  for 
space  diversity  has  already  been  discussed  in  the  previous  chapter  and  illustrated  in 
Figure  8.  It  is  assumed  that  perfect  CSI  is  known  at  the  receiver  end.  In  MRC,  the 
received  signal  is  multiplied  with  the  complex  conjugate  of  CSI  and  all  diversity 
receptions  are  added  to  form  the  decision  variable.  The  outputs  of  the  integrators  of  the 
MRC  for  each  diversity  reception  can  be  represented  as 

Yl=AThHlXk  +  Nl 

Yl  =  ATbH\Xk  +  N? 

Y1  -  Rb^k-^k  |  jyl'  (3.13) 

k  A  k 

Y2  =  EhHlXk  +  2, 

A 


where  Eh  is  the  average  energy  per  bit,  Nk  =  A  TbN\ ,  and  Nk  =  A  TbNk . 

After  the  integrator  stage,  the  random  variable  is  multiplied  by  the  complex 
conjugate  of  CSI  and  the  resulting  random  variable  can  be  represented  as 


EbHi(Hiyxk 

A 

A 


+  (HiTK 


+  (H2JNI 


EbH\:Xk 

A 


z;  = 


E„  H  X  X 


A 


(3.14) 


The  decision  variable  Zk  at  the  output  of  combiner  stage  can  be  expressed  as 

=  Z\  +  Z\ 

=  Et  (|if‘|2  +  \Hlf)  A  +  +  (Hi)*  N'2 ' 
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For  the  final  decision  of  the  estimated  received  data,  the  decision  variable  is 
compared  with  predetennined  threshold  levels.  For  the  case  of  BPSK,  the  real  part  of  Zk 
is  compared  with  a  threshold  level  of  zero  to  decide  each  received  binary  bit.  [2,  15] 

D.  SUMMARY 

This  chapter  discussed  the  designs  of  the  MDDM  transmitter  and  receiver  at  the 
block  level.  The  function  of  each  block  was  discussed  with  the  representation  of  signals 
before  and  after  each  block  operation.  During  the  discussion,  two  assumptions  were 
made,  i.e.,  perfect  knowledge  of  CSI  at  the  receiver  and  constant  channel  response  for  the 
duration  of  the  OFDM  symbol.  The  next  chapter  analyzes  the  MDDM  transmitter  and 
receiver  model  for  the  AWGN  channel  with  and  without  flat  fading. 
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IV.  ANALYSIS  AND  SIMULATION  OF  MULTICARRIER  DELAY 
DIVERSITY  MODULATION  SCHEME 


In  this  chapter,  the  performance  (bit  error  probability)  analysis  and  simulation  of 
MISO  and  MIMO  systems  with  multicarrier  delay  diversity  is  presented.  The  MDDM 
transmitter  and  receiver  as  described  in  Chapter  III  were  simulated  in  Matlab.  The  system 
is  analyzed  with  an  AWGN  channel  with  and  without  Rayleigh  fading.  The  analytical 
and  simulated  results  are  compared  with  a  SISO  BPSK  system.  To  establish  a  fair 
comparison  of  the  SISO  system  and  the  MIMO  system  with  MDDM,  the  power 
transmitted  for  both  systems  is  assumed  equal.  To  facilitate  better  understanding  and 
precise  representation  of  signals  from  transmit  to  receive  antenna,  this  thesis  will  adhere 
to  the  following  notation: 

•  Eb  represents  the  average  energy  per  bit  at  the  receiver  of  the  SISO 
system, 

•  Es  represents  the  average  energy  per  MIMO  BPSK  symbol  per  receiving 
antenna, 

•  Ts  represents  the  OFDM  symbol  duration, 

•  Ts  represents  the  MIMO  BPSK  symbol  duration  and  T  =  Tb 

•  £k  represents  the  real  part  of  decision  variable  Zk . 

Energy  per  bit  is  an  important  parameter  for  comparison  and  is  defined  as 

£.=J„'W \dt  =  AX  (4.1) 

where  k  is  the  time  index  of  the  BPSK  symbol  (before  serial  to  parallel  and  after 
parallel  to  serial  conversions)  and  Xk  is  the  BPSK  modulated  symbol  as  defined  in  Table 
2. 


Input  bit  at  time  k 

Output  Symbol  Xk 

0 

A 

1 

-A 

Table  2  BPSK  Modulation  Scheme  (After  Ref.  [8]) 
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For  simulation  and  analysis,  the  amplitude  A  =  1  is  maintained.  Equation  (4.1) 
shows  that  the  energy  per  bit  for  BPSK  modulation  is  the  same  for  bit  0  and  bit  1 .  Thus 
the  transmitted  power  for  the  SISO  system  can  be  written  as 

Pmo  =  ^  ^  (4.2) 

1b 

The  probability  of  bit  error  Pb  for  a  baseband  equivalent  SISO  BPSK  system  in 
discrete  domain  for  correlation  demodulator  can  be  represented  as  [15,  16] 


(4.3) 


The  bit  error  probability  of  the  SISO  system  as  given  by  Equation  (4.6)  will  serve  as  a 
benchmark  for  comparison  of  the  MISO  and  MIMO  systems  with  MDDM  in  the  AWGN 
channel. 

In  this  thesis  the  number  of  transmit  antennas  is  two  and  both  the  transmit 
antennas  transmit  equal  power.  The  total  energy  per  bit  transmitted  for  MISO  and  MIMO 
systems  using  BPSK  symbols  can  be  defined  as 

E\  =  PTS  (4.7) 
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and 


SISO 


-  p  -  p 

—  1  MISO  —  1  MIMO  ' 


(4.8) 


Therefore,  the  energy  transmitted  per  antenna  P  can  be  given  as 


P'  = 


p 

rSISO 

I. 


p 

rSISO 

2 


(4.9) 


Now,  Equation  (4.7)  is  rewritten  as 


,  ,  ,  P  T 

J7j  —  P  P  —  1  SISO±s 


E\=~- 

2 


SIMULATION  OF  MDDM  TRANSMITTER 


(4.10) 


The  block  diagram  of  the  MDDM  transmitter  model  is  shown  in  Figure  9.  The 
MDDM  transmitter  was  simulated  in  MATLAB  with  equivalent  baseband  BPSK  in  the 
discrete  time  domain.  The  simulation  was  implemented  with  one  sample  for  each  BPSK 
symbol.  The  MDDM  transmitter  scheme,  as  mentioned  in  Chapter  III  with  two 
transmitting  antennas,  was  simulated  without  added  guard  interval,  D/A  converter  and  RF 
modulator  blocks  to  facilitate  the  simulations.  Equal  power  was  transmitted  from  both  the 
antennas.  To  achieve  the  same  total  power  transmission  as  that  of  a  single  antenna  BPSK 
transmitter,  the  signal  at  each  branch  of  MIMO  transmitter  was  multiplied  with  a  gain 
factor  of  g  for  nonnalization.  The  transmitted  energy  per  symbol  for  BPSK  is  same 
whether  a  binary  1  or  0  is  transmitted.  With  the  gain  factor  g  the  energy  per  symbol  is 
represented  as 


g,=^\gXi\dt  =  g‘AX- 


(4.11) 


Substituting  Equation  (4.10)  into  Equation  (4.1 1)  yields 

iWL  =  (/ax  • 

2 


(4.12) 


Using  Equation  (4.2),  Equation  (4. 12)  can  be  rewritten  as 
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(4.13) 


9 2  = 


9  = 


2 

1 

V2' 


The  gain  factor  g  =  1  /  >/2  is  maintained  for  all  subsequent  simulations  to  transmit  the 
same  power  as  that  of  a  single  antenna  BPSK  system.  Thus  the  effective  amplitude  of  the 
BPSK  modulated  signal  for  each  antenna  is  given  as 


lyl  A=1 

1  fc|  V2 


B. 


SIMULATION  OF  MDDM  RECIVER 


(4.14) 


The  block  diagram  of  the  MDDM  receiver  model  is  shown  in  Figure  12.  The 
MDDM  receiver  was  also  simulated  in  MATLAB  with  equivalent  baseband  BPSK  in  the 
discrete  time  domain.  The  MDDM  receiver  simulation  was  implemented  as  discussed  in 
Chapter  III  except  for  the  RF  demodulator,  the  A/D  converter  and  the  Remove  Guard 
blocks.  The  receiver  was  configured  as  a  MISO  system  with  only  one  receive  antenna 
and  as  a  MIMO  system  with  two  receive  antennas  and  as  a  MIMO  system  with  three 
receive  antennas.  After  multicarrier  delay  diversity  demodulation,  the  space  diversity 
receptions  of  MIMO  systems  were  combined  by  using  the  optimum  MRC  technique  as 
discussed  in  Chapter  II.  After  the  BPSK  correlation  demodulator  as  illustrated  in  Figure 
14,  £k  (the  real  part  of  the  random  variable  Zk)  was  compared  with  the  threshold  level 
according  to  Table  3.  For  analysis  and  simulation  purposes  it  is  assumed  that  the 
transmitter  and  receiver  frequencies  are  synchronized  and,  in  the  case  of  MIMO  systems, 
all  the  diversity  receptions  are  also  synchronized. 
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Figure  14  BPSK  Correlation  Demodulator  for  MIMO  System  with  MDDM 


Output  Binary  Bit 

C->0 

0 

C*<  o 

1 

Table  3  Demodulation  of  BPSK  Signal  (After  Ref.  [8]). 

C.  SIMULATION  AND  PERFORMANCE  ANALYSIS  OF  MDDM  IN  AWGN 

The  perfonnance  of  MDDM  was  first  simulated  and  analyzed  in  AWGN  only.  In 
this  simulation  and  analysis  no  fading  is  assumed.  This  means  that  each  channel  response 
coefficient  hlj  equals  one.  A  block  diagram  of  the  MIMO  system  in  AWGN  with  three 
receiving  antennas  is  illustrated  in  Figure  15.  Both  the  MDDM  transmitter  and  receiver 
are  collapsed  into  one  block  each  to  facilitate  presentation.  In  this  simulation  two  gain 
blocks  each  with  normalizing  gain  factor  of  g  are  shown  at  the  outputs  of  the  transmitter 
for  each  transmitting  antenna.  The  AWGN  channel  blocks  add  white  Gaussian  noise  to 
the  signal  at  respective  receiving  antennas.  The  noise  power  is  increased  progressively 
with  each  simulation  run  to  calculate  the  bit  error  rate  at  different  signal  to  noise  ratios. 
Bit  error  rate  is  calculated  by  comparing  the  input  binary  data  stream  at  the  input  to  the 

transmitter  bk  and  output  binary  data  stream  at  the  output  of  the  receiver  bk  . 


53 


54 


where  N  is  the  FFT  size  or  the  number  of  subcarriers,  Nk  is  AWGN  in  the  frequency 
domain  and  (1  +  e~j2>ck/N)  is  the  effective  channel  response  for  MDDM  in  AWGN 
channel.  To  facilitate  the  presentation,  the  phase  angle  <j)k  is  defined  as 

(j)k  =  2nk  /  N.  (4.20) 

Multiplying  Equation  (4.19)  by  the  complex  conjugate  of  the  effective  channel  response, 
the  result  is  given  by 

Zk  =  (1  +  ejA )  (1  +  e+jA  )Xl  +  (1  +  e+jA  )Nk .  (4.2 1) 

Using  Euler’s  identity 

ej 6  =  cos(6*)  +  j  sin(0)  (4.22) 

Equation  (4.21)  is  rewritten  as 

+  =  (2  +  2co s(h))Xl  +  (1  +  (4.23) 

Equation  (4.23)  defines  a  complex  Gaussian  random  variable  Zk  due  to  AWGN. 
For  demodulation  of  BPSK  data  ,  £k  =  Re  j  Zk }  is  compared  with  the  threshold  as  given 
in  Table  3.  If  correlation  demodulator  conditions  are  assumed  then 


r  i  n  r  * 

E  [Zk  |  "  0"  was  transmitted}  =  E  <j  —  J  —j=  (2  +  2  cos  </>k )  +  Nk  (l  +  e+lA  } 


+  o  L 


dt 


A 


=  [2 +  2  cos  (&.)]-^r 
=  V2[l  +  cos(4)]A. 


(4.24) 


Since  this  expected  value  is  real,  it  follows  that 

^ k  =  E  |Re  [Zk  }  |  "0"  was  sent}  =  y/2  [l  +  cos(^t )]  A  (4.25) 


The  variance  of  £k  =  Re  { Zk }  is  only  due  to  the  variance  of  real  part  of  the  noise 
component  (1  +  e 1 1<k  )N k  [15,  16]  and  is  represented  by  cri 


55 


crj  =  E  (Re[(l  +  e*«)JVt])1  - (e{B*[(1  +  e+»)JVj 

(4-26) 

CT?  =E  (Re[(l  +  e*«R.])j 


where  E  {  }  represents  the  expected  value.  r/k  =  Re  [(1  +  e+3A  )Nk  ]  is  the  real  part  of  a 
complex  Gaussian  random  variable  with  zero  mean.  rjk  and  ///  can  be  written  as 

n„  =  Re[(l  +  e**)JVt]  =  ^((1  +  +  (1  +  e*)Nl) 

nl  =  (Re  [(1  +  e*»  )Nt  ]f  =  {  i  ((1  +  e*«  )Nt  +  (1  +  e*  )N‘t  )T 

V2  J  (4.27) 

U((l  +  e+3A  fN'l  +  (1  +  e-*  )27V;2  +  2(1  +  e-* )  (1  +  e+*  ) 

=  I((l  +  e^fNl  +  (1  +  e“*  W  +  4(1  +  cos(^))Ayv;). 

where  (  )*  represents  the  complex  conjugate.  Substituting  Equation  (4.27)  into  Equation 
(4.26)  gives 

of  =  eE((1  +  e“*f  (AT,)2  +  (1  +  e*f  (N;f  +  4(1  +  cos(A)WV;)} 

1  ’  (4.28) 

=  j((l  +  e*'*)2  E{(l\g2}  +  (1  +  e-’*)!  E{(JV;)2}  +  4(1  +  cosfo))E{Ar,lv;}). 

As  discussed  earlier,  Nk  is  a  complex  Gaussian  random  variable  with  zero  mean, 
therefore 

E{(JVt)2}  =  E{(JV;)2}  =  0.  (4.29) 

Substituting  Equation  (4.29)  into  Equation  (4.28)  gives 

ct!  =[l  +  cos(*)]E{ArtJV;}.  (4.30) 

Using  Equation  (2.29),  Nk  is  given  as 

N- 1 

Nk  =  (4.31) 

m= 0 

Substituting  Equation  (4.31)  into  Equation  (4.30)  yields 
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N- 1 


crj  =  [l  +  cos(4)]Ej^ 

[_m= 0 
N-l  N- 1 

=  [1  +  cos(^,)]£^E 

/=0  m= 0 
N-l  N-l 

=  [l  +  cos(4)]^^E 


N-l 


^-jZnkn/N  z  n*e+j2M/N 


\ninrn}e 


\ninm)e 


i=o  J 

-j2nkl  /  N  j2nkm  /  N 

-  j27ik{l-m)  /  N 


(4.32) 


/=0  m= 0 


Samples  of  AWGN  at  different  time  instants  are  IID  Gaussian  random  variables  with 
zero  mean.  Therefore,  they  are  uncorrelated  and  the  variance  can  be  represented  as 


E  {nmnp}  = 


^E{nm}E  jn*j  if  m*p 
E{|nm|2}  if  rn  =  p 

0  if  rn  p 

if  rn  =  p. 


(4.33) 


cr 


Substituting  Equation  (4.33)  into  Equation  (4.32)  gives 


J(^)]ZEKX, 


-  [l  +  COS 

=  [l  +  cos(^)]7Vcr^ 


(4.34) 


where  <72n  is  the  noise  power  of  the  OFDM  symbol  in  the  time  domain  before  the  FFT 

operation.  If  correlation  demodulator  conditions  are  assumed,  the  noise  power  for  OFDM 
symbol  durations  is  given  by 

J\T  00  ^ 

<  =  I  |tf(/)|  (4-35> 

-CO 

where  N0  /  2  (Watts  per  Hertz)  is  the  power  spectral  density  of  the  real  part  of  noise  and 
H(f)  is  the  frequency  response  of  the  integrator  [15,  16].  Using  Parsevel’s  theorem  [19], 
Equation  (4.35)  converts  to 
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(4.36) 


El 
2  Tt‘ 


Substituting  Equation  (4.36)  into  Equation  (4.34)  gives 

2  _  [l  +  cos(&)]MV0 
f  2  Ts 


(4.37) 


Substituting,  Ts  from  Equation  (2.33)  and  since  for  BPSK  T'  =  Tb ,  Equation  (4.37) 
yields 


.2  =  [l  +  cos(4)]##0  =  [l  +  cos(&)]iV0 
f  2  NTb  2  Tb 


(4.38) 


Using  Equations  (4.24)  and  (4.38),  the  bit  error  probability  Pb  conditioned  on  the  value 
(j)k  can  be  expressed  as 


Pb  =  Pr \^k  <  0  |  "0"  was  sent} 


=  Pr 


Ck  -  c 


r+ 

<-  _  hk  |  II  Q  III 


a 


=  Pr 


Ck  -  c 


<y  c 


=  Q 

=  Q 


bk 


=  Q 


y/2A(  l  +  cos(^)) 


^[l  +  cos  ($,.)]  No  /  2Th 
2A^Tb  [1  +  cos($ 


p 

f 

=  Q 

7 

V 

Eb  [1  +  cos(&) 


N, 


(4.39) 


The  number  of  data  symbols  per  OFDM  frame  is  192  as  discussed  in  Chapter  III. 
The  indices  of  these  data  symbols  are  given  as  k  =  9, 10, ...,100, 156, 157,..., 255  .  Thus, 
the  average  bit  error  probability  Pb  is  given  by 


58 


(4.40) 


The  simulation  was  conducted  with  increasing  Eb  /  N0  for  100,000  OFDM 
frames  or  19.2  million  data  symbols  per  Eb  /  No  value. 

The  simulated  bit  error  rate  (BER)  for  the  MISO  system  is  plotted  in  Figure  16 
where  L  =  2  is  the  number  of  transmit  antennas  and  J  =  1  is  the  number  of  receive 
antennas.  The  theoretical  probability  of  bit  error  as  obtained  in  Equation  (4.40)  is  also 
plotted.  For  comparison  of  performances,  this  figure  also  includes  the  theoretical 
probability  of  bit  error  for  a  baseband  equivalent  SISO  system.  The  simulated  results 
follow  the  theoretical  results  very  closely.  The  MISO  system  with  MDDM  performs 
better  than  the  SISO  system  for  lower  Eb  /  N0  values  and  the  performance  of  the  MISO 
system  is  poorer  than  that  of  SISO  system  for  Eb  /  No  greater  than  6.5  dB. 
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Figure  16  Results  of  MDDM  MISO  System  in  AWGN 
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2.  Performance  Analysis  of  MIMO  System  with  Two  Transmit  and  Two 
Receive  Antennas 

For  the  MDDM  MIMO  system  with  two  receive  antennas,  it  is  assumed  that  the 
signal  receptions  at  both  the  antennas  are  uncorrelated  and  signals  are  received  at  the 
same  time  without  any  relative  delay.  The  received  signals  at  both  receive  antennas  are 
given  by 


=  x„ 


=  x„ 


+  xl  +  n 
+  xl  +  n 


1 

m 

2 

m 


After  the  FFT  operation  the  signals  can  be  represented  as 

Rl=Xl+Xl+Ni 

K  =  A  +  x\  +  n\. 


Substituting  Equation  (2.31)  and  (4.19)  into  Equation  (4.42)  gives 

R\  =  X\  +  Xle~j2xk/N  +  Nl 
Rl  -  Xl  +  Xle-j2*k/N  +  N; 

A={l  +  e-*)Xl+Nl 
A=(l  +  e-*)X1k+N2i. 


(4.41) 


(4.42) 


(4.43) 


Following  the  derivation  of  Equation  (4.23),  Equation  (4.43)  converts  to 

2\  =  Rl(  1  +  e'* )  =  [2  +  2  cost*)]  +  (1  +  e**  )N\ 
Z\  =  7J|(1  +  e**)  =  [2  +  2cos(*)]x;  +  (1  +  e**)J S\. 


(4.44) 


Now  both  space  diversity  receptions  are  combined  to  form  random  variable  Zk  which  is 
given  by 


zk  =  z\  +  zl 

=  [4  +  4  cost*)]  Xl  +  (1  +  e*>*  )jv;  +  (1  +  )N; 

=  4  [1  +  cost*)]  +  (1  +  c*®  )K  +  (1  +  e*®  )Nl 
=  2j2A  [1  +  cos(^)]  +  (1  +  e+jA  )Nl  +  (1  +  e+jA  )N2k. 


(4.45) 


If  correlation  demodulator  conditions  are  assumed,  then 
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r  1  n  ] 

E{Zj,  |  "0"  was  transmitted}  =  E  j—  J  Zkdt  ^ 


=  2~j2A(l  +  cos  ($,.)). 


(4.46) 


The  variance  of  £k  =  Re  { Zk }  is  only  due  to  the  variance  of  real  part  of  noise 

components  [15,  16].  Noise  components  at  both  receive  antennas  are  IID  Gaussian 
random  variables.  Therefore  the  total  variance  of  the  sum  of  two  noise  components  is  the 
sum  of  their  individual  variances  [18].  The  total  noise  variance  is  given  as 

<  =  +  <4-47) 

where  C,'  =  He  \Z): }  and  ^  =  Re  \Z'k } .  Following  the  derivation  of  noise  variance 
derived  in  Equation  (4.37),  Equation  (4.47)  is  represented  as 

[1  +  cos (</>k)]N0  (  [l  +  cos(^)]i\^ 

(4.48) 


2  T 


2  T 


[l  +  cos(4)].¥0 


and  <j_  is  given  as 


^  = 


(1  +  cos  {</>k))No 


(4.49) 


Following  the  derivation  of  Equation  (4.39),  the  bit  error  probability  Pi  conditioned  on 
the  value  </>k  is  given  by 


h=q 

=  Q 

=  Q 


(  2^2  [l  +  cos(^.)]^ 


^[l  +  cos^J]  Na/Tb 


8  [l  +  cos(^.)]  A2Tb 


b  J 
\ 


(4.50) 


8[l  +  cos(^,)]E; 


N„ 
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Following  the  derivation  of  Equation  (4.40)  the  average  probability  of  bit  error  Pb  for  a 
MDDM  MIMO  system  with  two  transmit  and  two  receiving  antennas  is  given  by 


Pb  = 


192 


(  100 

f  f 

Z<2 

j 

\^k= 9 

1 8(1  +  cos  (^,  ))E), 
N 


1  255 

f  r 

+  Eo 

J- 

J  k= 156 

Iv 

1 8(l  +  cos(^))E) 
N 


(4.51) 


JJ 


The  simulation  was  conducted  for  the  same  number  of  OFDM  frames  or  data 
symbols  as  mentioned  in  the  previous  section.  The  simulated  bit  error  rate  (BER)  for  the 
MIMO  system  is  plotted  in  Figure  17.  The  simulated  results  follow  theoretical  results 
very  closely. 


Figure  17  Results  of  MDDM  MIMO  System  in  AWGN 


3.  Performance  Analysis  of  MIMO  System  with  Two  Transmit  and 
Three  Receive  Antennas 

Similarly,  in  this  section,  the  performance  analysis  of  MDDM  MIMO  system  with 
two  transmit  and  three  receive  antennas  is  discussed.  Exploiting  the  symmetry  of  the 
receiver  design,  the  derivation  of  results  is  mainly  based  on  the  results  presented  in  the 
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last  two  sections.  It  is  also  assumed  that  the  signal  receptions  at  all  antennas  are 
uncorrelated  and  signals  are  received  at  the  same  time  without  any  relative  delay.  The 
received  signals  at  the  three  receive  antennas  are  given  by 


112  1 

r  =  £  +  £  +  n 


2  1,2,2 
r,  =  £„,  +  £.„  +  ri 


1  1  2  2 
rl  =  xt,  +  xt  +  nt 


After  the  FFT  operation  the  signals  are  represented  as 

Rlk=Xl+Xl+Nl 
Rl=X\  +  Xl+Nl 
Rl=Xl  +  Xl+Nl 


(4.52) 


(4.53) 


Based  on  the  discussion  in  previous  sections,  the  components  of  the  variables  from  all  the 
space  diversity  receptions  are 

Zl  =  J?2(l  +  e**)  =  [2  +  2cos(<4)]X2  +  (1  +  e'*)N[ 

Zl  =  Rl(  1  +  )  =  [2  +  2  cos(  A)]  X,'  +  (1  +  )N;  (4.54) 

Zl  =  Rl(  1  +  e*»)  =  [2  +  2cos(A)]x;  +  (1  +  e**)Nl 

Now  all  the  space  diversity  receptions  are  combined  to  fonn  random  variable  Zk  which 
is  given  by 

zk  =  zl  +  zl  +  z3k 

(4.55) 

=  6[1  +  cos(A)]x;  +  (1  +  e*)Nl  +  (1  +  e*)Nl  +  (1  + 


If  correlation  demodulator  conditions  are  assumed,  then 


f  i  Tb 

E  [Zk  |  "  0"  was  transmitted}  =  E  l  —  J  ^2>/2v4  (1  +  cos  <f>k )  +  Nk  (l  +  e+:A  ) 

1^6  o 


dt 


(4.56) 


=  2\[2A\l  +  cos(^,)  . 


The  variance  of  £k  =  Re  { Zk }  is  only  due  to  the  variance  of  real  part  of  noise 
components  [15,  16].  Noise  components  at  all  the  receive  antennas  are  IID  Gaussian 


63 


random  variables.  Therefore  the  total  variance  of  the  sum  of  three  noise  components  is 
the  sum  of  their  individual  variances  [18].  The  total  noise  variance  is  given  as 


'Ck 


cl 


a 


Ck 


(4.57) 


where  ■£[.  =  Re { Z\  j ,  =  Re { Z\  j  and  =Re| ZkJ.  Following  the  derivation  of 

noise  variance  derived  in  Equation  (4.37),  Equation  (4.57)  is  represented  as 

3[l  +  cos  (</>k)\N0 


<  = 


2  T 


(4.58) 


and  <j:  is 


= 


|3[l  +  cos(4)]jVo 


2  T 


(4.59) 


Following  the  derivation  of  Equation  (4.39),  the  bit  error  probability  1\  conditioned  on 
the  value  (j>k  is  given  by 


H  =  Q 
=  Q 

=  Q 


3>/2  [l  +  COs(^.)]y4 


yj3[l  +  cos(<f>k)]N0/2Tb 
|l2[l  +  cos  {<f>k)\A2Tb 

V  -AT 


[12[l  +  cos(4)]Efc 
N„ 


(4.60) 


The  average  probability  of  bit  error  Pb  for  this  case  can  be  obtained  by 


P,  = 


192 


100 

f  [ 

ZQ 

V 

k=9 

V 

vV 

|12[1  +  cos(A  )]£'»''  2“  ^ 


N. 


k= 156 


12[l  +  cos(^)]Eft 

7V„ 


JJ 


(4.61) 


The  simulation  was  conducted  for  the  same  number  of  OFDM  frames  or  data 
symbols  as  mentioned  in  the  previous  section.  The  simulated  bit  error  rate  (BER)  for  the 
MIMO  system  is  plotted  in  Figure  18.  The  simulated  results  follow  theoretical  results 
very  closely. 
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A  comparison  of  Equations  (4.40),  (4.51)  and  (4.61)  clearly  reveals  that  the 
improvement  factor  increases  linearly  with  increasing  number  of  receiving  antenna.  Thus, 
the  average  probability  of  bit  error  in  AWGN  channel  for  MDDM  scheme  discussed  in 
this  thesis  with  J  receiving  antenna  can  be  represented  in  general  form  as 


Pb  = 


1 


192 


(  100 

(  r 

£q 

i 

\k=S 

Iv 

/4J(1  +  cos(0k))Eb 
N 


)  255 

(  r 

+ 

j 

J  k= 156 

Vi 

/4J(l  +  cos  {<f>k))Eb 
N 


(4.62) 
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D.  SIMULATION  AND  PERFORMANCE  ANALYSIS  OF  MDDM  IN  A 
MULTIPATH  FADING  CHANNEL 

Next,  the  MDDM  scheme  was  simulated  and  analyzed  with  the  effect  of  a  slow 
fading  frequency  nonselective  channel.  The  multipath  model  for  the  slow  Rayleigh  fading 
channel  was  discussed  in  Chapter  II.  A  multipath  fading  channel  is  inserted  in  Figure  15 
from  each  transmit  antenna  to  each  receive  antenna.  A  block  diagram  of  the  MIMO 


0  2 

Figure  18 
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system  with  slow  fading  frequency  nonselective  channel  plus  AWGN  is  illustrated  in 
Figure  19.  It  is  assumed  that  the  each  channel  frequency  response  remains  constant  for  an 
OFDM  symbol  duration.  Therefore  the  channel  is  slow  fading.  The  simulated  channel 
response  is  the  complex  sum  of  two  IID  Gaussian  random  sources.  The  block  model  of 
channel  response  is  illustrated  in  Figure  20.  In  simulation,  each  channel  response  is 
multiplied  with  the  transmitted  signal  from  the  respective  antenna  and  input  to  the 
AWGN  block  of  respective  receive  antenna. 


Figure  19  Simulation  of  MDDM  MIMO  System  in  Multipath 
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Figure  20  Simulation  of  Channel  Response 
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Each  Gaussian  random  source  generates  a  zero  mean  random  variable.  In  this 
simulation  model,  the  variance  of  each  Gaussian  source  is  equal  and  is  set  to  a  value  such 
that  the  variance  of  the  complex  sum  of  these  Gaussian  random  sources  is  one.  Using 
Equation  (2.52),  the  channel  response  from  transmit  antenna  l  to  receive  antenna  j  is 
represented  as 


hlj  =  \hlj\ej0,i 

=  X  +  jY 


(4.63) 


where  X  and  Y  are  IID  Gaussian  random  variables  with  zero  mean  and  equal  variance. 
The  variance  of  the  complex  sum  can  be  given  as 


a 


0 


=  E{hl3hli*}-{^E{hl3}^ 

=  E{(X  +  jY)(X  +  jY)*} 
=  E{X2  +Y2} 


9  2 

=  °x  +  °Y 


=  1. 


(4.64) 


Thus  the  variance  of  each  Gaussian  random  source  was  set  as 


(4  =<4=1/2.  (4.65) 

1.  Performance  Analysis  of  MISO  System  with  Two  Transmit  and  One 
Receive  Antenna 

At  the  receive  antenna,  the  signal  rm  is  written  as 

rm  =  hlxxxm  +  h21x2m  +  nm  .  (4.66) 

where  h 11  and  h 21  are  channel  responses  from  transmit  antenna  1  and  2  to  the  receive 
antenna  respectively.  These  channel  responses  are  assumed  independent.  After  the  FFT 
operation,  the  signal  is  written  as 

FFT[rJ  =  FFTlfcX  +  hnx2m  +  nj.  (4.67) 

All  the  channel  responses  are  considered  constant  for  the  duration  of  an  OFDM  symbol 
for  the  slow  fading  channel.  Now,  Equation  (4.67)  can  be  written  as 

Rk=hnXl+h21X2k+Nk.  (4.68) 
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Substituting  Equation  (2.3 1)  ,  (4.16)  and  (4.20)  into  Equation  (4.68)  yields 


Rk  =  hllX\  +  h21X\e*  +  Nk 
Rk=(hn+h21e-^)Xi+Nk 


(4.69) 


where  Nk  is  AWGN  in  the  frequency  domain  and  (//  '  +  hne  J*  )  is  the  effective  channel 

response  for  the  MDDM  in  the  multipath  channel.  In  Equation  (4.69)  h11  is  a  random 
variable  with  Rayleigh  distributed  amplitude  and  phase  uniformly  distributed  on  . 

Thus,  multiplying  h 21  by  e  1<!h  does  not  change  the  statistics  of  the  random  variable  and 
their  product  is  represented  as  h21  .  Now,  the  effective  channel  response  can  be  written  as 

H1  =  hu  +  h2V  (4.70) 


Using  the  MRC  receiver  as  shown  in  Figure  14,  Equation  (4.69)  is  multiplied  by 
the  complex  conjugate  of  the  effective  channel  response  and  the  and  random  variable  Zk 
is  represented  as 

Zk=H1H1*X1k+H1*Nk.  (4.71) 


Now,  Equation  (4.71)  can  be  written  as 

=  (, hn  +  h2V )(hn  +  h2V )*  X\  +  ( h 11  +  h2V )*  Nk  (4.72) 

=  (Kt  +  \hn'f  +  2(ft“Af  +  h'QX))xl  +  (A11  +  ft”')’  Nh 

where  h!J  =  hj  +  jhlQ  ,  i.e.  h\  and  fig  are  the  inphase  and  quadrature  components  of  the 
respective  channel  responses. 

For  a  fixed  set  of  channel  responses  Equation  (4.72)  represents  a  complex 
Gaussian  random  variable  Zk  due  to  AWGN  [16].  To  demodulate  the  BPSK  signal,  the 
real  part  of  the  decision  variable  Zk ,  £k  =  Re  { Zk } ,  is  compared  with  a  threshold  as 
given  in  Table  3.  If  correlation  demodulator  conditions  are  assumed,  then 
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Ck  =E{Re{Zk}  I  "0"  was  transmittedj 


b  0 
|2 


=  E  -J  \hn[  +\h21]  +2(h]X'  +KX)Xl+(hU  +h2V)  Nkdt\  (4.73) 


fc11|2+|^1,|2+2(W  +  W) 


V2 


The  variance  of  £k  is  only  due  to  the  variance  of  real  part  of  the  noise  component 


Hv>Nk  [15,  16]  and  is  represented  as 


{H 

f] 

|(Re| 

>x] 

)1 

(4.74) 


where  E  {  }  represents  the  expected  value.  rjk  is  the  real  part  of  noise  component  with 
zero  mean  which  can  be  written  as 


%  =  Re[ff‘X]  =  +  H'K)  (4.75) 

And 

nl  =  (Refff'-JV,])2  =  +  H'Nljj 

=  \{(HV)2  Nl  +  (H'f  (Nl)2  +  2/71/f1*JVX)  (4.76) 

=  i((/71*)Wf  +(H'f(Nlf  +2\H'\‘ NtNl). 


Substituting  Equation  (4.74)  into  Equation  (4.75)  yields 

<4  =E{i((ff,’)!jv‘2  +  X)2  Wf  +  2M1jvx)} 

=  +  (H'f  E{W)2}  +  2|7f‘f  EjJVX}). 

Nk  is  a  complex  Gaussian  random  variable  with  zero  mean,  therefore 

E{(Kf}  =  E{(Nlfr0. 

Substituting  Equation  (4.78)  into  Equation  (4.77)  yields 


(4.77) 


(4.78) 
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.2 


cr 


(k 


E{jvtAr;} 

f  + 1*21'!2  +2(W  +  W))E{^Jv;}. 


(4.79) 


Now,  Equation  (4.79)  can  be  written  as 
2  1 


= 


rM^T +2(k,x'+w) 

||A,1|s+|AJ11!+2(ft"Af  +W)] 


271 


^  = 


Now,  the  conditional  bit  error  probability  can  be  represented  as 
Pb  =  Pr{^  <  0  |  "0"  was  sent} 

_  Ck  ~  ^k  ^  _  Ck  I  I!  Q  I!  I 


471 


cr 
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cr 


(k 


_  pr  Cfc  Cl:  ^  Cfc  I  I!  Q  I!  I 


cr 


<T* 


cr 
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=  Q 


=  Q 


bk 


+\h’r\'  +2(h?h?'  +KK))A/Ji 

)(|A''f+|A”'f+2(/!;'Af  +W))NJ4Tl 


=  Q 


27l2Tt  ( 

h11 

2 

+ 

h21' 

2  +  2(W  +  >W)) 

V  N0 

(4.80) 


(4.81) 


or 


pM  =  Q 


2  EbPi 


N. 


(4.82) 


where  /?,  is  given  as 

A=|ft“|!  +  |A2,f  +2(AX'+W).  (4.83) 

Equation  (4.82)  depicts  the  bit  error  probability  of  a  MDDM  MISO  system  conditioned 
on  the  random  variable  Pl .  The  perfonnance  of  an  MDDM  MISO  system  over  a 
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frequency  nonselective,  slowly  fading  Rayleigh  channel  with  MRC  can  be  obtained  by 
taking  the  average  of  the  conditional  bit  error  probability  over  all  possible  values  of  the 
random  variable  /?, .  Now,  the  unconditional  bit  error  probability  can  be  given  as 


0 

00 

=  \Q{sf2E^jW„)t,i(Pl)dpi 

0 


(4.84) 


where  /B  ( /?, )  is  the  probability  density  function  for  /?, . 

The  random  variable  /?,  as  defined  in  Equation  (4.83)  is  a  sum  of  two  central  chi 

squared  random  variables  of  degree  two,  h  and  h  ,  and  two  products  of  Gaussian 

random  variables,  2  hfh21'  and  2  h^h^' .  These  central  chi  squared  random  variables  and 

Gaussian  random  variables  are  also  correlated.  The  probability  density  function  of  the 
random  variable  /?,  and  integral  of  Equation  (4.84)  were  computed  numerically. 

The  random  variable  Pl  was  generated  in  Matlab  by  taking  one  million  samples 

each  of  channel  response  h 11  and  h2V  .  The  effective  channel  response  was  multiplied  by 
its  complex  conjugate.  The  distribution  of  /?,  over  its  range  was  obtained  by  evaluating  a 
histogram  of  500  equally  spaced  bins  with  the  help  of  the  Matlab  function  hist .  The 
distribution  of  the  data  was  interpolated  by  cubic  spline  data  interpolation  with  the  help 
of  the  Matlab  function  spline .  After  data  interpolation,  this  data  distribution  curve  was 
nonnalized  to  make  the  total  area  under  the  distribution  curve  equal  to  one.  Then,  the 
probability  over  all  possible  values  was  interpolated  by  using  the  Matlab  function  spline  . 
At  the  end,  the  average  probability  of  bit  error  was  computed  numerically  over  all  the 
possible  values  of  the  random  variable  Pl .  The  Matlab  code  is  presented  in  Appendix  B. 
[21] 

The  simulated  and  the  numerically  computed  theoretical  results  of  a  MISO  system 
over  frequency  nonselective  slowly  fading  Rayleigh  channel  are  shown  in  Figure  21.  For 
the  comparison  of  the  perfonnances,  the  theoretical  bit  error  rate  (BER)  of  the  BPSK 
SISO  system  as  obtained  in  Equation  (2.127)  is  also  plotted  along  with  the  bit  error 
probability  obtained  from  the  simulation.  The  simulated  results  are  close  to  the 
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theoretical  results  for  lower  values  of  Eb  /  No  and  for  higher  values  of  Eb  /  No  the 
theoretical  results  are  more  optimistic.  This  may  be  due  to  the  use  of  the  Matlab  function 
spline  to  interpolate  the  probability  distribution  function  of  the  random  variable  . 
Further  analysis  could  be  carried  out  to  determine  the  causes  of  this  deviation.  It  was  not 
investigated  further  in  this  work.  Figure  21  clearly  shows  that  the  performance  of  the 
MDDM  MISO  system  is  better  than  the  SISO  system. 


Figure  2 1  Results  of  MDDM  MISO  System  in  Slow  Rayleigh  Fading  Channel 


2.  Performance  Analysis  of  MIMO  System  with  Two  Transmit  and  Two 
Receive  Antennas 

In  this  section,  the  performance  analysis  of  a  MIMO  system  with  two  transmit  and 
two  receive  antennas  is  presented.  The  simulation  model  and  analysis  are  mainly  based 
on  the  previous  section.  At  each  of  the  receive  antennas,  the  received  signals  are  given  by 


1  ill  1  .  7  21  'A  . 

r„=h  xm  +  h  +  n. 


21  2 


2  7,12  1  .  7,22  2  .  2 

r„,  =  h  x,„  +  lx  x.„  +  n 


(4.85) 
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where  hlj  represents  channel  response  from  transmit  antenna  /  to  receive  antenna  j  and 
r;'  is  the  received  signal  at  receive  antenna  j .  It  is  considered  that  all  the  channel 
responses  are  independent  and  the  reception  of  the  signals  at  both  the  antennas  is 
uncorrelated.  After  the  FFT  operation  the  signals  are  represented  as 


FFT[r-]  =  FFT[/,X+kX+<] 

(4.86) 

FFT  [r^]  =  FFTpX  +  hVm  +n2m] . 

All  the  channel  responses  are  considered  constant  for  the  duration  of  an  OFDM  symbol. 
Now,  Equation  (4.86)  can  be  written  as 


Rt  =  hnX\  +  h21X2k  +  Nl 
Rl  =  h12Xl  +  h22X2k  +  N2k. 


(4.87) 


Substituting  Equations  (2.31)  ,  (4.16)  and  (4.20)  into  Equations  (4.87)  yields 


Rl  =  huXl  +  IrX'c  *  +  Nl 
R2k  =  h12Xlk  +  If  2  X  *  c  4  +  N2k 

(4.88) 

Rl  ={hu  +h21e-jA)Xl  +  Nl 
Rl  =  (h12  +  h22e~jA)Xl  +  N2k 

where  NJk  is  AWGN  in  the  frequency  domain  at  receive  antenna  j ,  (h  "  +  h2'e~:i^  j  is  the 
effective  channel  response  at  receive  antenna  1  and  (VJ  +  h22e~j^  j  is  the  effective 
channel  response  at  receive  antenna  2.  In  Equation  (4.88),  all  hlj  (l  =  1,2  and  j  =  1,2) 
are  random  variables  with  Rayleigh  distributed  amplitude  and  phase  uniformly 
distributed  on  (-#,#].  Thus,  multiplying  h 21  and  h22  by  e  J'A  does  not  change  the 
statistics  of  both  the  random  variables  and  their  respective  products  are  represented  as 
h2V  =  h2 1  e“^':  and  h22'  =  h22e  jVA  .  Now,  the  effective  channel  responses  are  written  as 


H1  =  hu  +  h21' 
H2  =  h 12  +  h22'. 


(4.89) 
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Using  the  MRC  receiver  as  shown  in  Figure  14,  Equation  (4.88)  is  multiplied  by 
the  complex  conjugate  of  the  effective  channel  response  and  the  random  variables  for 
respective  diversity  receptions  are  given  by 

Z\  =  I!'!!1*  Xj„  +  H'*Nlk 
z\  =  H2H2*X\  +  H2*Nl 

l  |  ]|2  i  v  ,  (4-90) 

z\  =  | Hx\  Xl  +  H  Nl 

Zl  =  Iff2 12  Xl  +  H2*N2k. 


Now  Equations  (4.90)  can  be  expressed  as 


Z\  =  (|A“f  +  |ft”f  +  2 {hfhf  +  h^))xl  +  (ft11  +  h“'j  Nl 
Z\  =  (|A'f  +  |A“'f  +  2 (hfhf  +  h^hf  ))xl  +  (A12  +  A22')’  JV,2 


(4.91) 


where  h\  =  Re|/r'j  and  =Im|h,J|  are  the  inphase  and  quadrature  components  of 
the  respective  channel  responses. 

After  combining  both  the  space  diversity  receptions  in  the  MRC  receiver  the 
decision  variable  Zk  is  given  as 


zk  =  Zl  +  Zl 

=  | H1 12  Xl  +  | H2  \  Xl  +  H^Nl  +  H2*N2k  (4.92) 

=  (|Af‘|!  +  |if2|2)x;  +  ff‘’jv;+Af2*iv2- 


For  a  fixed  set  of  channel  responses  Equation  (4.92)  represents  a  complex 
Gaussian  random  variable  Zk  due  to  AWGN  [16].  To  demodulate  the  BPSK  signal,  the 

real  part  of  the  decision  variable  Zk ,  C,k  =  Re  j  Z,  j  ,  is  compared  with  a  threshold  as  given 
in  Table  3.  If  correlation  demodulator  conditions  are  assumed,  then 
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c  =E{R*{Z1}|"0"  was  transmitted} 


=  E  — 


kfW 


+  \m~\Xl  +  H1*N\  +  H2*Nl 


\dt 


_A_ 

V2  L 

M 

V2 


h 11  f  +  |/t21'  f  +  2  (hfhf'  +  hfhf )  +  |  h12  f  +  |/i22’  f  +  2  (hfhf  +  h22h22’ ) 


(4.93) 


where 


ft  =  |/i“|3  +  |A”f  +  |A“f  +  |Ai2'|2  +  2(h)X'  +  K'h" 

2(W+A“Af)- 

The  variance  of  £k  is  only  due  to  the  variance  of  real  part  of  the  noise  components 
H1*^  and  H2*Nl  [15,  16].  Noise  components  at  both  the  receive  antennas  N\  and  Nf. 
are  IID  complex  Gaussian  random  variables.  Thus,  the  variance  of  their  sum  is  the  sum  of 
their  individual  variances  and  is  given  as 

<=d;+<4'  (4-95> 

Following  the  results  derived  in  the  previous  section  from  Equation  (4.74)  to  (4.80),  cri 
can  be  expressed  as 


(4.94) 


2 


=  \ | H1  f  E {NlN? }  +  \\H2 12  E {NlN2; } 


_  N0 


47], 


+ 


|A11f+|A2"|-+2(A“Af  +AX) 

|A‘f +|A!1f+2(A“Af  +A“Af) 


Substituting  Equation  (4.94)  into  Equation  (4.96)  yeilds 


2  _  -^oA 

4T 


= 


M 

1  47], 


(4.96) 


(4.97) 


Now,  substituting  Equation  (4.93)  and  (4.97)  into  Equation  (4.3),  the  conditional  bit  error 
probability  is  represented  as 
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(4.98) 


f 

PM  =  Q 


v 

/ 

=  Q  . 

V 


M/j 

^N0/ATbJ 


N. 


Equation  (4.98)  represents  the  bit  error  probability  of  an  MDDM  MIMO  system 
with  two  transmit  antennas  and  two  receive  antennas  conditioned  on  a  random  variable 
ft 2 .  The  performance  of  the  MDDM  MIMO  system  over  a  frequency  nonselective  slowly 
fading  Rayleigh  channel  with  MRC  can  be  obtained  by  taking  the  average  of  the 
conditional  bit  error  probability  over  all  possible  values  of  random  variable  /?2 .  Now  the 
unconditional  bit  error  probability  can  be  given  as 

oo 

(4.99) 

=  J  Q^fiJ  N„)  4(A)  dfi, 

0 

where  /B  ( /f )  is  the  probability  density  function  for  /?2  . 

The  random  variable  /?2  as  defined  in  Equation  (4.94)  is  a  sum  of  central  chi 
squared  random  variables  of  degree  two  and  products  of  Gaussian  random  variables.  The 
central  chi  squared  random  variables  and  Gaussian  random  variables  are  also  correlated. 
The  probability  density  function  of  the  random  variable  /f  and  integral  of  Equation 
(4.99)  were  computed  numerically,  using  the  same  algorithm  as  discussed  in  previous 
section. 

The  simulated  and  the  numerically  computed  theoretical  results  of  the  MIMO 
system  with  two  transmit  antennas  and  two  receive  antennas  over  a  frequency 
nonselective  slowly  fading  Rayleigh  channel  are  shown  in  Figure  22.  For  the  comparison 
of  the  performances,  this  figure  also  includes  the  theoretical  and  the  simulated  bit  error 
probability  of  a  BPSK  SISO  system.  The  simulation  results,  for  low  Eb  /  N0  values, 
deviate  slightly  from  the  theoretical  results.  The  theoretical  results  are  more  optimistic  for 
low  Eb  /  No .  The  reason  for  this  slight  deviation  may  be  the  same  as  discussed  in  the 
previous  section,  however,  this  was  not  investigated  in  this  work. 
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3.  Performance  Analysis  of  MIMO  System  with  Two  Transmit  and 
Three  Receive  Antennas 

In  this  section  the  performance  analysis  of  a  MIMO  system  with  two  transmit  and 
three  receive  antennas  is  presented.  The  scheme  for  building  the  simulation  model  and 
the  analysis  is  the  same  as  discussed  in  last  two  sections.  At  the  receive  antennas,  the 
received  signals  are  given  by 


it  =  hLLxl  +  h2lxl  +  n. 


2  7  12  1  ,  7  22  2  ,  2 

r  =  h  xm  +  h  +  n 


rl=h13xi+h23xi  +nl 


After  the  FFT  operation  the  signals  are  written  as 


(4.100) 
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(4.101) 


FFT[^,]  =  FFTpX  +  Hnxl  +  n'„] 

FFT[^]  =  FFT[/iX  +  AX  +  n;] 

FFT[r*]  =  FFT[A“4  +  Kax‘„  +  <]. 

All  the  channel  responses  are  considered  fixed  for  the  duration  of  an  OFDM  symbol. 
Now,  Equation  (4.101)  can  be  written  as 

R\  =  huX\  +  K2lXl  +  N\ 

Rl  =  h12Xl  +  h22  Xl  +  N'l  (4. 1 02) 

Rf.  =  hl3X\  +  h23X2k  +  TV3. 

Substituting  Equations  (2.31),  (4.16)  and  (4.20)  into  Equations  (4.102)  yields 

R\  =  hnX\  +  h2lX\e-^  +  N\ 

R;  =  hl2X\  +  h22Xle~jA  +  N2k 
R!  =  hl3X\  +  h23Xle-jA  +  TV3 

(4.103) 

R\  =  (hu  +  h21e-jA)Xl  +  Nl 
R;  =  (h21  +  h22e-jA)Xl  +N2k 
Rl  =(h13  +h23e^)Xl  +  Nl 

where  Nk  is  AWGN  in  the  frequency  domain.  (//"  +  h21e  M  ),  (/r  +  h22e  7'A  )  and 
(h13  +  h23e~j^)  are  the  effective  channel  responses  at  receive  antennas  1,  2  and  3 
respectively.  In  Equation  (4.103),  all  hl 3  (l  =  1,2  and  j  =  1,2,3)  are  random  variables 
with  Rayleigh  distributed  amplitude  and  phase  uniformly  distributed  in  Thus, 

multiplying  h21 ,  h22  and  h23  by  e  does  not  change  the  statistics  of  these  random 
variables  and  their  respective  products  are  represented  as  h2V  =  h21e  J"/'i  ,  h22'  =  h22e~3^ 
and  h23  =  h23e  7'A  .  Now,  the  effective  channel  responses  are  given  as 

H 1  =  hu  +  h2V 

H2=hv2  +  h22'  (4.104) 

H3  =  h 13  +  h23'. 
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Using  MRC  receiver  as  shown  in  Figure  14,  Equation  (4.88)  is  multiplied  by  the 
complex  conjugate  of  the  effective  channel  response  and  the  random  variable  Z[  at 
respective  diversity  receptions  are  written  as 


Z\  =  H'H'*X'k  +  H'*N\ 

Z\  =  H2H2*X\  +  H2*N'l 
Z\  =  H3H3*Xl  +  H3*N3k 

z\  =  \H3\2  Xl  +  h1*n\ 

Z2k  =  \H2(  Xl  +  H2*N2k 
Z\  =  \H3f  xl  +  H3*Nl 


(4.105) 


Now,  Equations  (4.105)  can  be  represented  as 


Z\ 

K 

zl 


|ft”f + H”'f +2  (h)X‘ + KK')]xl  +  0“ + '>"')■ 

\h'^ +\h^ +2(hfh?  +AJV  )]x;  +0“ 

|am|!  +  |A"f  +2(  hfhf  +  W)]x;  +  (ha  +  h*)'  Ni 


(4.106) 


where  h!]  =  Re|/;/'|  and  hg  =Im|h,J|  are  the  inphase  and  quadrature  components  of 
the  respective  channel  responses. 

After  combining  all  the  space  diversity  receptions  in  the  MRC  receiver  the 
decision  variable  Zk  is  given  as 


zk=zl  +  z2k+z3k 

=  | H1 12  Xl  +  | H2 12  Xl  +  | H3 \2  Xl  +  HX*N\  +  H2*N'l  +  H3*N3k  (4. 107) 

=  (I#1)2  +  \H2f  +  |i73|2)  Xl  +  Hx*Nl  +  H2*N2k  +  H3*Nl 


For  a  fixed  set  of  channel  responses  Equation  (4.107)  represents  a  complex 
Gaussian  random  variable  Zk  due  to  AWGN  [16].  To  demodulate  the  BPSK  signal,  the 
real  part  of  the  decision  variable  Zk ,  £k  =  Re  { Zk } ,  is  compared  with  a  threshold  as 
given  in  Table  3.  If  correlation  demodulator  conditions  are  assumed,  then 
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C=E{Re{Zt}|"0"  was  transmitted} 

=  e|  — J|"(|iy1|2  +  |tf2|2  +  \H3\2)jXl  +  Hl*N\  +  H2*N2h  +  H3*N3k  dt 

=  {\H1f  +\H2f  +\H3\2)JXl  (4.108) 

A  H/t11!2  +|/t21'|2  +2(W  +  h1Q%v)  +  \h12\2  +\h22'\2  + 

2  (hfhf  +  h%hf )  +  |/t13  f  +  \h 23’  |2  +  2  ( W  +  W ) 

V2 

where 

yS3  =  |/i“f  +  +  |ft"f  +  +  |A1S|3  +  |/12S'|2  +2  (hi'hf  +  h^X' 

2(fc“Af+W)+2(fc,X'  +  W)- 

The  variance  of  is  only  due  to  the  variance  of  real  part  of  the  noise  components 
Hl*N\,  H2*N;  and  II2'  N 2  [15,  16].  Noise  components  at  all  the  receiving  antennas 
Nl ,  Nl  and  N2  are  IID  complex  Gaussian  random  variables.  Thus,  the  variance  of  their 
sum  is  the  sum  of  their  individual  variances  and  is  represented  as 

d.=d;+d;+<-  <4'110> 

Following  the  results  derived  in  previous  section  from  Equation  (4.74)  to  (4.80),  od  can 
be  written  as 


(4.109) 


°l  = 


E  {NIK }  +  K  E{JV|JV|*}  +  \H’f  E  {jV’JVf }) 

N0  \hU\  +  |^21 1  +  2(hfhf  +  IiqIiq  )  +  |/t12|  +  \h22  |  + 

4^  2  (hfhf  +  hX')  +  |/r13|2  +  l^'f  +  2(  W  +  W) 


(4.111) 


Substituting  Equation  (4. 109)  into  Equation  (4.111)  yields 

2  _  NqP3 


<  = 


= 


(4.112) 
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Now,  substituting  Equation  (4.108)  and  (4.112)  into  Equation  (4.3),  the  conditional  bit 
error  probability  is  written  as 


PM)  =  Q 


'  APJ-J2  ' 


»  J 


=  Q 


\2EbP3 


N 


O  J 


(4.113) 


Equation  (4.113)  shows  the  bit  error  probability  of  the  MDDM  MIMO  system 
with  two  transmit  antennas  and  three  receive  antennas  conditioned  on  the  random 
variable  J3:i .  The  performance  of  the  MDDM  MIMO  system  over  a  frequency 
nonselective,  slowly  fading  Rayleigh  channel  with  MRC  can  be  obtained  by  taking  the 
average  of  the  conditional  bit  error  probability  over  all  possible  values  of  random 
variable  fl,  .  Now,  unconditional  bit  error  probability  can  be  represented  as 

oo 

P>  =  \PM)h,(h)M 

(4.114) 

o 

where  /B  ( /?, )  is  the  probability  density  function  for  //, . 

The  random  variable  /?.,  as  defined  in  Equation  (4.109)  is  also  a  sum  of  central 
chi  squared  random  variables  of  degree  two  and  products  of  Gaussian  random  variables. 
The  central  chi  squared  random  variables  and  the  Gaussian  random  variables  are  also 
correlated.  The  probability  density  function  of  the  random  variable  /?.,  and  integral  of 
Equation  (4.99)  were  computed  numerically,  using  the  same  algorithm  as  discussed  in 
previous  section. 

The  simulated  and  the  numerically  computed  theoretical  results  of  the  MIMO 
system  (with  2  transmit  antennas  and  3  receive  antennas)  over  a  frequency  nonselective 
slowly  fading  Rayleigh  channel  are  shown  in  Figure  23.  For  comparison  of  the 
performances,  this  figure  also  includes  the  theoretical  and  the  simulated  bit  error 
probability  for  a  BPSK  SISO  system.  This  figure  shows  that  for  low  Eb  /  N0  values  the 
theoretical  results  are  more  optimistic  and  for  high  Eb  /  N0  values,  the  simulation  results 
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are  more  optimistic.  This  deviation  of  simulation  results  from  the  theoretical  results  may 
also  be  due  to  the  reasons  as  discussed  in  previous  sections.  Further  analysis  can  be 
carried  out  to  investigate  the  exact  causes  of  this  trend  as  a  future  work;  this  was  not 
investigated  further  in  this  work. 


Figure  23  Results  of  MDDM  MIMO  System  in  Slow  Rayleigh  Fading  Channel 

E.  SUMMARY 

In  this  chapter,  the  MDDM  MISO  and  MIMO  systems  as  discussed  in  Chapter  III 
were  examined  with  the  AWGN  channel  with  and  without  slow  Rayleigh  fading.  These 
systems  were  simulated  in  Matlab  and  the  theoretical  results  were  shown  in  comparison 
with  the  simulated  results.  The  simulated  results  follow  the  theoretical  results  closely.  In 
the  next  chapter  the  results  of  the  thesis  are  summarized  and  areas  for  follow  on  work  are 
presented. 
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V.  CONCLUSION 


The  goal  of  this  thesis  was  to  design  and  analyze  a  MIMO  system  with 
multicarrier  delay  diversity  modulation.  MDDM  was  incorporated  with  orthogonal 
frequency  division  multiplexing  and  cyclic  delay  diversity.  A  MISO  system  and  MIMO 
systems  were  designed  using  the  MDDM.  They  were  simulated  in  Matlab  with  BPSK  and 
were  tested  in  a  progressive  manner,  first  in  an  AWGN  channel  and  then  in  a  multipath 
fading  channel  with  AWGN.  The  BER  performance  of  these  systems  was  analyzed  and 
compared  with  the  performance  of  the  simulated  systems.  The  simulated  performance 
results  and  theoretical  analysis  results  were  compared  with  the  conventional  SISO  system 
performance. 

A.  RESULTS 

The  simulated  MDDM  MISO  and  MIMO  systems  achieved  BER  perfonnance 
results  consistent  with  the  theoretical  analysis  in  an  AWGN  channel  with  and  without 
multipath  fading.  The  perfonnance  metric  of  bit  error  probability  versus  Eb  j  N0  (energy 
per  bit  to  noise  power  spectral  density  ratio)  was  used.  To  establish  a  fair  comparison,  the 
transmitted  power  and  data  rate  for  SISO,  MISO  and  MIMO  systems  were  equal.  The 
comparison  of  perfonnances  in  AWGN  showed  that  the  MDDM  MISO  system 
performed  better  than  the  SISO  system  for  low  Eb  /  N0  with  up  to  a  6.5  dB  performance 
gain  and  performed  worse  for  higher  Eb  /  No  values.  The  performances  of  the  MDDM 
MIMO  systems  were  better  than  that  of  SISO  system  for  all  values  of  Eb  / N0.  MIMO 
systems  with  two  receive  antennas  and  three  receive  antennas,  outperformed  a  SISO 
system  by  1.7  dB  and  3.4  dB  performance  gain.  In  a  multipath  fading  channel  with 
AWGN,  MISO  and  MIMO  systems  were  able  to  achieve  significant  advantage  over  a 
SISO  system.  The  improvement  in  performances  of  the  MDDM  MISO  and  MIMO 
systems  over  the  SISO  system  can  be  exploited  for  military  and  civilian  applications. 

B.  RECOMMENDATION  FOR  FUTURE  RESEARCH 

There  are  four  areas  identified  for  future  research.  First,  the  MDDM  MIMO 
system  performance  can  be  analyzed  in  different  fading  conditions  as  described  by 
another  fading  model  in  place  of  the  slow  frequency  nonselective  Rayleigh  fading 
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channel.  Good  choices  for  alternate  fading  models  include  slow  Ricean  frequency 
nonselective  fading,  slow  Nakagami  frequency  nonselective  fading,  and  various  forms  of 
frequency  selective  fading. 

In  addition  to  the  analysis  performed  in  this  thesis,  higher  order  modulation 
schemes  in  lieu  of  BPSK,  i.e.,  M  PSK  and  M  QAM,  can  be  incorporated  and  analyzed  to 
evaluate  the  performance  of  the  MDDM  MIMO  systems. 

Next,  forward  error  correction  (FEC)  coding  with  interleaving  can  be 
incorporated  into  the  MDDM  MIMO  system  with  a  maximum  likelihood  sequence 
estimation  receiver.  The  performance  can  be  analyzed  theoretically  and  compared  with 
simulated  performance. 

Lastly,  the  MIMO  system  analysis  herein  assumes  perfect  CSI  at  the  receiver.  If 
this  infonnation  about  the  channel  can  be  extracted  at  the  receiver,  then  it  can  be  sent  to 
the  transmitter.  The  transmitter  can  use  CSI  to  more  efficiently  transmit  power  to  the 
receiver  through  channels  that  are  less  faded  than  others.  Techniques  that  control  the 
transmitter  can  be  investigated  and  perfonnance  of  the  system  analyzed. 
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APPENDIX 


MATLAB  CODES 


********* 


%  Number  of  subcarriers 
%  Data  symbol  in  an  OFDM  frame 
%  Total  number  of  data  for  the  simulation 
%  high  frequency  carriers 
%  low  frequency  carriers 
%  All  pilot  sybmols  set  to  ones 
%  DC  component  not  transmitted 


This  appendix  includes  the  Matlab  code  for  simulating  the  MDDM  scheme  in  an 
AWGN  channel  with  and  without  frequency  nonselective  Rayleigh  fading.  This  also 
includes  the  theoretically  computed  results  for  the  BER  performance  of  the  MDDM. 

A.  SIMULATION  OF  THE  MDDM  IN  AN  AWGN  CHANNEL 

%********  MDDM  Simulation  with  AWGN  only 

clear  all 

num_frame  =  100000; 
fftsize  =  256; 
data_size  =192; 
t_data  =  data_size*num_frame; 
hfc  =  27; 

Ifc  =  28; 
pilot=ones(1 ,8); 
dc_comp  =  0; 
for  n=1:fftsize 

p_shift(n,1)=exp(i*2*pi*(n-1)/fftsize);  %  Phase  shift  due  to  cyclic  delay 
end  %  afer  fft  operation 

EbNo_dB=0:12; 
for  nn=1  :length(EbNo_dB) 

SNR=EbNo_dB(nn)+  3;  %  converting  SNR  to  EbNo 
t_err_SISO=0;  t_err_2T1R=0;  t_err_2T2R=0;  t_err_2T3R=0; 

BER=0;  totalerr=0; 

%**********  Simulating  input  data  and  BPSK  modulation 
for  mm  =1:num_frame 
r  =  rand(1,data_size); 
for  m=1:data_size 
if  r(1,m)<  0.5 

info(1,m)=0;  symbol(1,m)=1; 
else 

info(1,m)=1;  symbol(1  ,m)=-1; 
end 


************* 


end 


Creating  OFDM  Block 


************ 


Pre_IFFT_blk=[  dc_comp  pilot  symbol(1:92)  zeros(1,(hfc+lfc))  symbol(93:data_size)]; 
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0/o*************  simulating  Trnasmitted  signal************* 
txsigl  =  ifft((Pre_IFFT_blk).');  %  taking  IFFT  after  S/P  coversion 
txsig2  =  circshift(txsig1 ,1 );  %  adding  Cyclic  Delay  Diversity 

Tx_signal  =  1/sqrt(2)*(txsig1+txsig2); 

%*********  Simulating  Receivs  ********* 

Rx_sig  =  awgn(symbol,SNR,  measured');  %  Received  signal  for  SISO  BPSK 
Rx_sig_1  =  awgn(Tx_signal, SNR, 'measured');  %  Received  signal  by  Ant.  1 
Rx_sig_2  =  awgn(Tx_signal, SNR, 'measured');  %  Received  signal  by  Ant.  2 
Rx_sig_3  =  awgn(Tx_signal, SNR, 'measured');  %  Received  signal  by  Ant.  3 
Rx_2T  1 R  =  fft(Rx_sig_1 );  %  2  Tx  1  Rx  Antennas 

Rx_2T2R  =  fft(Rx_sig_1+Rx_sig_2);  %  2  Tx  2  Rx  Antennas 
Rx_2T3R  =  fft(Rx_sig_1+Rx_sig_2+Rx_sig_3);  %  2  Tx  3  Rx  Antennas 
Rx_data_2T  1 R  =(Rx_2T1  R  *(1  +p_shift)); 

Rx_data_2T2R  =(Rx_2T2R.*(1+p_shift)); 

Rx_data_2T3R  =(Rx_2T3R.*(1  +p_shift)); 

Rx_data_2T  1 R  =(Rx_data_2T  1 R).';  %  converting  P/S 

Rx_data_2T2R  =(Rx_data_2T2R).'; 

Rx_data_2T3R  =(Rx_data_2T3R).'; 

%  Extracing  received  data  from  ofdm  block 

Rx_info_2T  1 R  =  [Rx_data_2T1R(10:101)  Rx_data_2T1R(157:fftsize)]; 
Rx_info_2T2R  =  [Rx_data_2T2R(10:101)  Rx_data_2T2R(157:fftsize)]; 
Rx_info_2T3R  =  [Rx_data_2T3R(10:101)  Rx_data_2T3R(157:fftsize)]; 

0^*******************  00pp|Q(jy |g^jpjg  BP3K  SiQDSl  ****************** 


for  k=1:data_size 

if  real(Rx_sig(1,k))  >  0 
R_info(1,k)=0; 
else  R_info(1,k)=1; 
end 

if  real(Rx_info_2T  1  R(1  ,k))  >  0 
R_info_2T1R(1,k)=0; 
else  R_info_2T  1  R(1  ,k)=1 ; 
end 

if  real(Rx_info_2T2R(1,k))  >  0 
R_info_2T2R(1  ,k)=0; 
else  R_info_2T2R(1  ,k)=1; 
end 

if  real(Rx_info_2T3R(1,k))  >  0 


%  demodulating  SISO  info 


%  demodulating  2T1R  info 


%  demodulating  2T2R  info 


%  demodulating  2T3R  info 
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R_info_2T3R(1  ,k)=0; 
else  R_info_2T3R(1  ,k)=1; 
end 
end 

%************  Probability  of  bit  error  Pb  ************* 

err_SISO  =  length(find(info-R_info));  %  calculating  errors 

err_2T  1 R  =  length(find(info-R_info_2T1R)); 

err_2T2R  =  length(find(info-R_info_2T2R)); 

err_2T3R  =  length(find(info-R_info_2T3R)); 

t_err_SISO  =  t_err_SISO+err_SISO;  %  Calculating  total  errors 

t_err_2T  1 R  =  t_err_2T  1  R+err_2T  1 R; 

t_err_2T2R  =  t_err_2T2R+err_2T2R; 

t_err_2T3R  =  t_err_2T3R+err_2T3R; 

end 

Pb_SISO(nn)  =  t_err_SISO/t_data;  %calculating  BER 
Pb_2T  1  R(nn)  =  t_err_2T1  R/t_data; 

Pb_2T2R(nn)  =  t_err_2T2R/t_data; 

Pb_2T3R(nn)  =  t_err_2T3R/t_data; 
end 

save  MDDM_AWGN_R  EbNo  dB  Pb_SISO  Pb_2T1R  Pb_2T2R  Pb_2T3R 

load  Theo_AWGN_R 

figure 

semilogy(EbNo_dB,Pb_SISO_T,'kl,EbNo_dB,  Pb_2T1  R_T,'r',EbNo_dB,  Pb_2T2R_T,... 
'b',EbNo_dB,  Pb_2T3R_T,,g,,EbNo_dB,Pb_SISO,'ko,,  EbNo_dB,Pb_2T1R,'rv',... 
EbNo_dB,Pb_2T2R,'b*',EbNo_dB,Pb_2T3R,,gd') 
grid  on 

xlabel('E_b/N_o  in  dB') 
ylabel('P_b') 

legend('SISO_T','2T  1  R_T','2T2R_T','2T3R_T',... 

'SISO_S','2T  1  R_S','2T2R_S','2T3R_S') 
title('Comparison  of  Theoretical  and  Simulated  Resulsts') 

B.  COMPUTING  THEORETICAL  BER  OF  THE  MDDM  IN  AN  AWGN 
CHANNEL 

%*******  Theoratical  BER  in  AWGN  Channel  ******** 

clear  all 

phase  =  linspace(0,2*pi*255/256,256); 
phi  =  [phase(10:101)  phase(157:256)]'; 
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comp  =  (1+cos(phi)); 

EbNo_dB=0:12; 

EbNo=10.A(EbNo_dB/10); 
for  n  =  1  :length(EbNo) 
temp  =  EbNo(n); 
for  k=1:length(phi) 

Q2T1R(k,n)  =  0.5*erfc(sqrt(2*temp*comp(k))); 

Q2T2R(k,n)  =  0.5*erfc(sqrt(4*temp*comp(k))); 

Q2T3R(k,n)  =  0.5*erfc(sqrt(6*temp*comp(k))); 
end 
end 

%*****  Theoretical  BPSK  ****** 

Pb_SISO_T=0.5*erfc(sqrt(EbNo)); 

Pb_2T  1  R_T  =  1/length(phi)*sum(Q2T  1 R); 

Pb_2T2R_T  =  1/length(phi)*sum(Q2T2R); 

Pb_2T3R_T  =  1/length(phi)*sum(Q2T3R); 

save  Theo_AWGN_R  EbNo  dB  Pb  SISO  T  Pb_2T1R_T  Pb_2T2R_T  Pb_2T3R_T 

load  M  D D M_AW G N_R 

figure 

semilogy(EbNo_dB,Pb_SISO_T,'k',EbNo_dB,  Pb_2T1  R_T,'r',EbNo_dB,... 
Pb_2T2R_T,'b',EbNo_dB,  Pb_2T3R_T,  g',EbNo_dB,Pb_SISO,'ko',... 

EbNo_dB,Pb_2T1R,'rv',  EbNo_dB,Pb_2T2R,'b*',EbNo_dB,Pb_2T3R,,gd') 
grid  on 

xlabel('E_b/N_o  in  dB') 
ylabel('P_b') 

legend(,SISO_T';2T1R_T,,,2T2R_T',,2T3R_T','SISO_S','2T1R_S','2T2R_S','2T3R_S') 
title('Copmarison  of  Theoretical  and  Simulated  Resulsts') 

C.  SIMULATION  OF  THE  MDDM  IN  FREQUENCY  NONSELECTIVE 
SLOW  FADING  RAYLEIGH  CHANNEL 

%*  Simulation  of  MDDM  MISO  and  MIMO  Systems  in  Frequency  Nonselective  Slow 
%  Fading  Rayleigh  Channel 
%  RAYLEIGH  FADING  CHANNEL 
clear  all 

num_frame=  100000; 

fftsize=256;  %  Number  of  subcarriers 

data_size=192;  %  Data  symbol  in  one  frame 

t_data  =  data_size*num_frame;  %  Total  number  of  data  for  the  simulation 
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%  high  frequency  carriers 
%  low  frequency  carriers 


%  All  pilot  sybmols  set  to  ones 
%  DC  component  not  transmitted 


hfc  =  27; 

Ifc  =  28; 
pilot=ones(1 ,8); 
dc_comp  =  0; 
for  n=1:fftsize 

p_shift(n,1)=exp(-i*2*pi*(n-1)/fftsize);  %  Phase  shift  due  to  cyclic  delay 

end  %  afer  fft  operation 

EbNo_dB=0:10; 

for  nn=1  :length(EbNo_dB) 

SNR=EbNo_dB(nn)+3; 

t_err_SISO_F=0;  t_err_2T1R_F=0;  t_err_2T2R_F=0;  t_err_2T3R_F=0; 

BER=0;  totalerr=0; 

%*********  Simulating  input  data  and  BPSK  modulation 
for  mm  =1:num_frame 
r=rand(1  ,data_size); 
for  m=1:data_size 
if  r(1,m)<  0.5 

info(1,m)=0;  symbol(1,m)=1; 
else 

info(1,m)=1;  symbol(1,m)=-1; 

end 
end 


********* 


^^  ************** ******  OFDM  Block  ************ 

Pre_IFFT_blk=[  dc_comp  pilot  symbol(1 :92)  zeros(1,(hfc+lfc))... 
symbol(93:data_size)]; 

0/o*************  simulating  Trnasmit  signal************* 

txsig1=  1/sqrt(2)*ifft((Pre_IFFT_blk).');  %  taking  IFFT  after  S/P  coversion 
txsig2=circshift(txsig1 ,1); 

%*********  Simulating  Received  s  i  g  n  a  I  %  %  %  %  %  %  %  %  %  %  %  %  %  %  % 

%  Frequency  Responses  of  Channels 

h=1/sqrt(2)*(randn(1  ,data_size)+j*randn(1  ,data_size));  %  SISO  system  channel 
h_1 1  =  1/sqrt(2)*(randn(1)+j*randn(1));  %  Tx  Ant.  1  to  RxAnt.  1 
h_21  =  1/sqrt(2)*(randn(1)+j*randn(1));  %Freq.  Response  of  Chan.  T.A.2  to  R.A.1 
h_12=1/sqrt(2)*(randn(1)+j*randn(1));  %Freq.  Response  of  Chan.  T.A.1  to  R.A.2 
h_22=1/sqrt(2)*(randn(1)+j*randn(1));  %Freq.  Response  of  Chan.  T.A.2  to  R.A.2 
h_13=1/sqrt(2)*(randn(1)+j*randn(1));  %Freq.  Response  of  Chan.  T.A.1  to  R.A.3 
h_23=1/sqrt(2)*(randn(1)+j*randn(1));  %Freq.  Response  of  Chan.  T.A.2  to  R.A.3 
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H_1  =h_1 1+h_21*p_shift;  %  Composite  channel  response  at  R.A.  1 
H_2=h_12+h_22*p_shift;  %  Composite  channel  response  at  R.A.  2 
H_3=h_13+h_23*p_shift;  %  Composite  channel  response  at  R.A.  3 
%********  Simulating  Recevers  ******* 

Rx_sig  =  awgn(symbol.*h, SNR, 'measured');  %  Received  signal  for  SISO  BPSK 
Rx_sig_1  =awgn((txsig  1  *h_1 1  +txsig2*h_21 ), SNR, 'measured'); 

%  simulating  Ant.#  1  received  signal  with  rayleigh  flat  fading  channel 
Rx_sig_2=awgn((txsig1*h_12+txsig2*h_22), SNR, 'measured'); 

%  simulating  Ant.#  2  received  signal  with  rayleigh  flat  fading  channel 
Rx_sig_3=awgn((txsig1*h_13+txsig2*h_23), SNR, 'measured'); 

%  simulating  Ant.#  2  received  signal  with  rayleigh  flat  fading  channel 

fft _ sig _ 1  =fft(Rx_sig_1 ); 

fft _ s  i  g_2=fft(  Rx_s  i  g_2 ) ; 

fft _ s  i  g_3=fft(  Rx_s  i  g_3 ) ; 

%Simulating  MRC  Receiver 

Vk_SISO=Rx_sig.*conj(h);  %SISO  system  reception  with  CSI 
Vk_1  =fft_sig_1  .*conj(H_1);  %Antenna  1  diversity  reception 
Vk_2=fft_sig_2.*conj(H_2);  %Antenna  2  diversity  reception 
Vk_3=fft_sig_3.*conj(H_3);  %Antenna  3  diversity  reception 
%  Combining  all  space  diversity  receptions 

Zk_2T  1  R=real((Vk_1 ).');  %Decision  variable  of  MISO  (2T 1 R)  and  P/S  conversion 

Zk_2T2R=real((Vk_1+Vk_2).');  %Decision  variable  of  MIMO  (2T2R)  and  P/S  conversion 
Zk_2T3R=real((Vk_1+Vk_2+Vk_3).');  %Decision  variable  of  MIMO  (2T3R)  and  P/S  conversion 
%  Extracing  real  part  of  decision  variable  from  ofdm  block 
Rx_info_2T  1  R_F  =  [Zk_2T1  R(1 0: 1 01 )  Zk_2T1  R(157:fftsize)]; 

Rx_info_2T2R_F  =  [Zk_2T2R(10:101)  Zk_2T2R(157:fftsize)]; 

RxJnfo_2T3R_F  =  [Zk_2T3R(10:101)  Zk_2T3R(157:fftsize)]; 

RxJnfo_SISO_F=real(Vk_SISO);  %BPSK  SISO  System  received  data 
^*******************  Demodulating  BPSK  signal  ****************** 
for  k=1:data_size 

if  Rx_info_SISO_F(1  ,k)  >  0  %  demodulating  SISO  info 

R_info_SISO_F(1  ,k)=0; 
else  R_info_SISO_F(1  ,k)=1 ; 
end 

if  Rx_info_2T  1  R_F(1  ,k)  >  0  %  demodulating  2T 1 R  info 

R_info_2T1R_F(1,k)=0; 
else  R_info_2T1R_F(1,k)=1; 
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%  demodulating  2T2R  info 


end 

if  Rx_info_2T2R_F(1  ,k)  >  0 
R_info_2T2R_F(1,k)=0; 
else  R_info_2T2R_F(1  ,k)=1 ; 
end 

if  Rx_info_2T3R_F(1  ,k)  >  0  %  demodulating  2T3R  info 

R_info_2T3R_F(1,k)=0; 
else  R_info_2T3R_F(1,k)=1; 
end 
end 

%%************  Probability  of  bit  error  Pb  ************* 

err_SISO_F  =  length(find(info-RJnfo_SISO_F));  %  calculating  errors 

err_2T  1  R_F  =  length(find(info-RJnfo_2T1R_F)); 

err_2T2R_F  =  length(find(info-RJnfo_2T2R_F)); 

err_2T3R_F  =  length(find(info-RJnfo_2T3R_F)); 

terrSISOF  =  t_err_SISO_F+err_SISO_F;  %  Calculating  total  errors 

t_err_2T1  R_F  =  t_err_2T1  R_F+err_2T1  R_F; 

t_err_2T2R_F  =  t_err_2T2R_F+err_2T2R_F; 

t_err_2T3R_F  =  t_err_2T3R_F+err_2T3R_F; 

end 

Pb_SISO_F(nn)  =  t_err_SISO_F/t_data;  %calculating  BER 
Pb_2T1R_F(nn)  =  t_err_2T1R_F/t_data; 

Pb_2T2R_F(nn)  =  t_err_2T2R_F/t_data; 

Pb_2T3R_F(nn)  =  t_err_2T3R_F/t_data; 
end 

save  MDDM_SRF_R  EbNo  dB  Pb_SISO_F  Pb_2T1R_F  Pb_2T2R_F  Pb_2T3R_F 

load  Theo_SRF_R 

figure 

semilogy(EbNo_dB,Pb_SISO_F_T,'k',EbNo_dB,  Pb_2T1R_F_T,'r',EbNo_dB,... 
Pb_2T2R_F_T,'b',EbNo_dB,  Pb_2T3R_F_T,'g',EbNo_dB,Pb_SISO_F,'ko',... 
EbNo_dB,Pb_2T1R_F,'rv',  EbNo_dB,Pb_2T2R_F,'b*',EbNo_dB,Pb_2T3R_F,,gd') 
grid  on 

xlabel('E_b/N_o  in  dB') 
ylabel('Pb') 

legend('SISO_T','2T  1  R_T',,2T2R_T,,'2T3R_T',,SISO_S','2T  1  R_S','2T2R_S','2T3R_SI) 
title('Copmarison  of  Theoretical  and  Simulated  Resulsts') 
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D.  COMPUTING  THEORETICAL  BER  OF  THE  MDDM  IN  FREQUENCY 
NONSLECECTIVE  SLOW  FADING  RAYLEIGH  CHANNEL 

%  Calculating  Theoretical  Results  By  Numerical  Methods 
global  WXN 
EbNo_dB  =  0:10; 

EbNo  =  10.A(EbNo_dB/10); 

%  creating  independent  channel  responses 

h_1 1  =  1/sqrt(2)*(randn(1 , 1 0A6)+j*randn(1 , 1 0A6));  %  Tx  Ant.  1  to  Rx  Ant.1 
h_21  =  1/sqrt(2)*(randn(1,10A6)+j*randn(1,10A6));  %  T.A.2  to  R.A.1 
h_1 2=  1/sqrt(2)*(randn(1 , 1 0A6)+j*randn(1 , 1 0A6));  %  T.A.1  to  R.A.2 
h_22=1/sqrt(2)*(randn(1 , 1 0A6)+j*randn(1 , 1 0A6));  %  T.A.2  to  R.A.2 
h_1 3=  1/sqrt(2)*(randn(1 , 1 0A6)+j*randn(1 , 1 0A6));  %  T.A.1  to  R.A.3 
h_23=1/sqrt(2)*(randn(1 , 1 0A6)+j*randn(1 , 1 0A6));  %  T.A.2  to  R.A.3 
H_1=h_11+h_21;  %*p_shift;  %  Effective  hannel  response  at  R.A.  1 
H_2=h_12+h_22;  %*p_shift;  %  Effective  channel  response  at  R.A.  2 
H_3=h_13+h_23;  %*p_shift;  %  Effective  channel  response  at  R.A.  3 
%  creating  random  variables  Beta 
Beta_1  =((H_1  ).*conj(H_1 )); 

Beta_2=((H_1  ).*conj(H_1  ))+((H_2).*conj(H_2)); 
Beta_3=((H_1).*conj(H_1))+((H_2).*conj(H_2))+((H_3).*conj(H_3)); 

[N 1  X1]=hist(Beta_1 ,500); 

[N2  X2]=hist(Beta_2,500); 

[N3  X3]=hist(Beta_3,500); 

%  Theoretical  BER  foR  SISO 
x  =  sqrt(2*EbNo./(1+2*EbNo)); 

Pb_SISO_F_T  =  0.5*(1-x); 

%  Estimating  PDF  Beta_1  and  Calculting  Theoretical  BER  for  2T1R 
X=X1;  N=N1 ; 

W=quad(@Func_Est_Cur,min(X1),max(X1));  %  Normalizing  factor  for  PDF 
Pr_Beta_1  =  quad(@Func_pdf_Beta,min(X1),  max(XI)) 
dBeta_1  =(min(X1  ):0.0001  :max(X1 )); 
for  i=1  :length(EbNo) 

Q_2Beta_1  =  0.5*erfc(sqrt(2*EbNo(i)*dBeta_1)); 

Pb_2T  1  R_F_T(i)  =  sum(Q_2Beta_1.*abs(Func_pdf_Beta(dBeta_1)))*0.0001 ; 
end 

%  Estimating  PDF  Beta_2  and  Calculting  Theoretical  BER  for  2T2R 
X=X2;  N=N2; 
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W=quad(@Func_Est_Cur,min(X2),max(X2));  %  Normalizing  factor  for  PDF 
Pr_Beta_2  =  quad(@Func_pdf_Beta,min(X2),  max(X2)) 
dBeta_2=(0:0.0001  :max(X2)); 
for  i=1  :length(EbNo) 

Q_2Beta_2  =  0.5*erfc(sqrt(2*EbNo(i)*dBeta_2)); 

Pb_2T2R_F_T(i)  =  sum(Q_2Beta_2.*abs(Func_pdf_Beta(dBeta_2)))*0.0001 ; 
end 

%  Estimating  PDF  Beta_3  and  Calculting  Theoretical  BER  for  2T3R 
X=X3;  N=N3; 

W=quad(@Func_Est_Cur,min(X3),max(X3));  %  Normalizing  factor  for  PDF 
Pr_Beta_3  =  quad(@Func_pdf_Beta,min(X3),  max(X3)) 
dBeta_3=(0:0.0001  :max(X3)); 
for  i=1  :length(EbNo) 

Q_2Beta_3  =  0.5*erfc(sqrt(2*EbNo(i)*dBeta_3)); 

Pb_2T3R_F_T(i)  =  sum(Q_2Beta_3.*abs(Func_pdf_Beta(dBeta_3)))*0.0001 ; 
end 

save  Theo_SRF_R  EbNo  dB  Pb_SISO_F_T  Pb_2T1R_F_T  Pb_2T2R_F_T  Pb_2T3R_F_T 

load  MDDM_SRF_R 

figure 

semilogy(EbNo_dB,Pb_SISO_F_T,'k',EbNo_dB,  Pb_2T1R_F_T,'r',EbNo_dB,  ... 
Pb_2T2R_F_T,'b',EbNo_dB,  Pb_2T3R_F_T,'g',EbNo_dB,Pb_SISO_F,'ko',... 
EbNo_dB,Pb_2T1R_F,,rv',  EbNo_dB,Pb_2T2R_F,'b*,,EbNo_dB,Pb_2T3R_F,'gd') 
grid  on 

xlabel('E_b/N_o  in  dB') 
ylabel('P_b') 

legend('SISO_T','2T  1  R_T','2T2R_T','2T3R_T','SISO_S','2T  1  R_S','2T2R_S','2T3R_S') 
title('Comparison  of  Theoretical  and  Simulated  Resulsts') 

E.  FUNCTIONS  TO  INTERPOLATE  PROBABAILITY  DISTRIBUTION 
FUNCTIONS 

function  y  =  Func_Est_Cur(x); 

%Title  :  Estimation  of  Data  Distribution  of  Random  Variable  Beeta 
%Author  :  Muhammad  Shahid,  Naval  Post  Graduate  School,  Septermber  2005 

% - 

%  y  =  Func_Est_Cur(x) 

% - 

%lnput  :  Random  variabel  Beeta 

%Output  :  Estimated  Data  Distribution  of  Beeta 
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% - 

global  WXN 
y=spline(X,N,x); 

function  PDF  =  Func_pdf_Beta(b_2); 

%Title  :  Intercalation  of  PDF  of  Random  Variable  Beeta 

%Author  :  Muhammad  Shahid,  Naval  Post  Graduate  School,  Septermber  2005 

% . 

%  PDF  =  Fun c_pdf_B eta ( b_2 ) ; 
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